ASYMPTOTIC FEYNMAN-KAC FORMULAE FOR LARGE 
SYMMETRISED SYSTEMS OF RANDOM WALKS 



By Stefan Adams 1 and Tony Dorlas 2 



(10 October, 2006) 



Abstract: We study large deviations principles for TV" random processes on the lattice 
Z rf with finite time horizon [0, f3] under a symmetrised measure where all initial and 
terminal points are uniformly given by a random permutation. That is, given a permu- 
tation a of N elements and a vector (xi, . . . , xn) of N initial points we let the random 
processes terminate in the points (x a n) , . . . , aww)) and then sum over all possible per- 
mutations and initial points, weighted with an initial distribution. There is a two-level 
random mechanism and we prove two-level large deviations principles for the mean of 
empirical path measures, for the mean of paths and for the mean of occupation local 
times under this symmetrised measure. The symmetrised measure cannot be written 
as any product of single random process distributions. We show a couple of important 
applications of these results in quantum statistical mechanics using the Feynman-Kac 
formulae representing traces of certain trace class operators. In particular we prove 
a non-commutative Varadhan Lemma for quantum spin systems with Bose-Einstein 
statistics and mean field interactions. 

A special case of our large deviations principle for the mean of occupation local 
times of N simple random walks has the Donsker- Varadhan rate function as the rate 
function for the limit N — > oo but for finite time (3. We give an interpretation in 
quantum statistical mechanics for this surprising result. 
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1. Introduction 

Let N random processes in continuous time on the lattice 7L d with initial distribution m G ty(Z d ), 
where is the set of probability measures on Z d , be given. We fix the time horizon as [0,(3]. In 

this paper we study large deviations for different functionals of the N random processes for large N 
under the symmetrised distribution 

1 N 

p S? = m E E ■•• E (!■!) 

Here ©at is the set of all permutations of N elements and the measure P^ j)X(r(j) is defined for any 
a G ©at and Xj G Z d , 1 < i < iV, as the conditional probability measure for the i-th random process 
starting at X{ with terminal location x a u\. That is, in (1.1) we have two mechanisms. First we 
draw uniformly a permutation and after that we pick N initial points which are permuted according 
to the chosen permutation to obtain N terminal points. Then these N initial and terminal points 
determine the N random processes. Finally we average over all permutations and initial points which 
are weighted with the given initial distribution m. Hence we prove two-level large deviations principles 
for the symmetrised distribution Pj^f. The symmetrised measure P^f itself is of interest because 
of the following reasons. 

The symmetrisation in (1.1) is described by the set of N pairs (x±, . . . ,x^; a^m, . . . , x a t N \) for any 
permutation a G &n and any x G 7j d . The mixing procedure for the second entry in these pairs 
has been studied both in [DZ92] and [Tr02], which were motivated from asymptotic questions about 
exchangeable vectors of random variables. [DZ92] studies large deviations for the empirical measures 
~k Si=i &Yi, where Y\, . . . ,Yn have distribution J Q fi(d6) Pffl for some distribution fi on some compact 
space ©, and the empirical measures are assumed to satisfy a large deviation principle under Pff for 
each 6. In [Tr02], a similar problem is studied: given a sequence of random vectors (Y-[ N \ . . . ,Y^) 
such that the empirical measures J2iLi ^ Y (N) sa ^ s fy a l &r g e deviation principle, another principle 

i 

is established for the process of empirical measures Y1a=\ $ x (n), where 

i 

fyW Y<> N A — J \ (vW v( N ) \ 

yy\ 1 , . . . , u\ N ) — ^ ^ \ cr(l)' ' ' ' ' cr(A^) / 
' ctGSjv 

Our second main motivation for studying the symmetrised distribution Pj^f stems from the applica- 
tion of Feynman-Kac formulae to express thermodynamic functions in quantum statistical mechanics. 
These thermodynamic functions are given as traces over exponentials of the Hamilton operator de- 
scribing the quantum system. There exist two kinds of elementary particles in nature, the Fermions 
and the Bosons. The state of a system of N Bosons is described by a symmetrisation procedure like in 
(1.1), whereas the state for Fermions is given with the corresponding anti-symmetrisation procedure. 
Thus one is lead to employ large deviation technique to study the large iV-limit for expectations with 
respect to the symmetrised distribution. We apply our main large deviation results in Section 2.2 and 
Section 2.3 to systems of Bosons in quantum statistical mechanics. 

We derive large deviations principles under the symmetrised distribution P^ ? for the empirical path 
measure Ljv, for the mean path Yn and for the mean of occupation measures Zn, all defined as functions 
of the N random paths . . . , £ (JV) : [0, f3) -► Z d , which are elements of the space Dp = D([0, f3\ 1 Z d ) of 
all functions uj : [0, 0\ — ► Z d , which are right continuous with left limits. The empirical path measures 



Ln = 7vE^« 
i=i 



(1.2) 
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are random elements in the set ty(Dp) of probability measures on Dp, the mean path 



1 N 

i=i 



(1.3) 



is a random element in Dp([0, 0\;M. d ) whereas the mean of the normalised occupation local times 



1 N 



i=i 



t(i) 

■f3 ' 



(1.4) 



is a random probability measure on Z d , where the normalised occupation local times are defined as 



which represent the relative time the i-th random process spends up to time (3 in the state z. 

Our large-deviation rate functions for the three principles are explicit in terms of variational prob- 
lems involving an entropy term (describing the large deviations of the permutations) and a certain 
Legendre transform (describing the large deviations of Ljy, Yn and Zjy, respectively, for a fixed permu- 
tation). These two parts in the variational formula for the rate function are due to the two- level large 
deviations, which has something in common with the multilevel large deviations studied in [DG94]. 
We draw a number of conclusions about variants of the principles, laws of large numbers and asymp- 
totic independence. Let us remark that all our large deviations results maybe obtained for random 
processes (Markovian or not) on any connected graph with finite or enumerable vertices. 

A first application of our large deviation results is given in Subsection 2.2, where we use the Feynman- 
Kac formula to represent the trace of any trace class operator restricted to the symmetric subspace 
of the iV-th tensor product of n x n complex matrix algebras as an expectation with a measure for N 
Markov processes on the index set {1, . . . m}. The trace class operator here is given by the Boltzmann 
factor e~@ h for any self-adjoint matrix h representing in quantum mechanics the Hamilton operator for 
a system of N quantum spins (lattice systems) for the inverse temperature (3, i.e. the time horizon of 
our random processes is given by the inverse temperature. In particular we derive the thermodynamic 
limit of the free energy, which is the trace of the Boltzmann factor, for a general class of mean-field 
interactions and thus we get a non-commutative version of Varadhan's Lemma with Bose-Einstein 
statistics, i.e. where the trace is restricted to the symmetric subspace. Here, an analysis of the 
variational formula for the rate function is achieved with L 2 techniques for the mean paths, which are 
embedding in the corresponding L 2 space. Non-commutative Varadhan's Lemmas have been studied 
in [CLR88] and [PRV89]. In [Dor96] a non-commutative central limit theorem under Bose-Einstein 
statistics has been proved for the case n = 2. Hence our results complement and extend these results. 

If we consider N simple random walks on Z^ conditioned to stay within a finite set A and replace 
the initial probability distribution m by the counting measure in A we are able to show that the 
corresponding rate function for the large deviations principle for the mean of occupation local times of 
the N random walks under the symmetrised measure (2.35) is given by the well-known Donsker- 

Varadhan rate function. The latter governs the large deviations principle for the occupation local time 
of a single random walk but for the limit (3 — > oo. This remarkable result has an interpretation for the 
cycle-structure given by our symmetrisation procedure, i.e. the appearance of cycles whose lengths 
grow like some potential of N. This long cycles are considered to be an order parameter for the 
occurrence of the Bose-Einstein condensation (BEC), a quantum phase transition solely driven by the 
symmetrisation procedure ([BCMP05],[DMP05]). Details of this interpretation are given at the end 
of Subsection 2.3. 




i = 




(1.5) 
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We consider this as a first step towards a rigorous understanding of large Boson systems at positive 
temperature /3, because the time horizon (3 represents the inverse temperature for the Feynman- 
Kac formulae. Future work will be devoted to the mutually interacting case. Interacting Brownian 
motions in trap potentials so far have been analysed without symmetrisation, in particular, systems 
for vanishing temperature in [ABK06a] and large systems of interacting motions for fixed positive 
temperature in [ABK06b]. In [AK06] some results for Brownian motions under the symmetrised 
distribution are obtained, which are not so general and which cannot be applied to mean field models. 

Let us make some remarks on related literature. We found a most interesting old work [Sch31] 
by Schrodinger, which is related to the pair probability method we applied in our large deviations 
principle. In [Sch31] Schrodinger raised the question of the most probable behaviour of a large system 
of diffusion particles in thermal equilibrium. Follmer [F688] gave a mathematical formulation of these 
ideas in terms of large deviations. He applied Sanov's theorem to obtain a large deviations principle 
for Ltv when B {1 \ B {2 \ . . . are i.i.d. Brownian motions with initial distribution m and no condition 
at time j3. The rate function is the relative entropy with respect to J Rd m(dx) F x o B^ 1 , where the 
motions start in x under P x . Then Schrodinger's question amounts to identifying the minimiser of 
that rate function under given fixed independent initial and terminal distributions. It turns out that 
the unique minimiser is of the form f Rd f Rd dxdy f(x)g(y) Fx, y ° B~ l , i.e. a Brownian bridge with 
independent initial and terminal distributions. The probability densities / and g are characterised 
by a pair of dual variational equations, originally appearing in [Sch31] for the special case that both 
given initial and terminal measures are the Lebesgue measure. 

An important work combining combinatorics and large deviations for symmetrised measures is 
[T6th90] . Toth [T6th90] considers N continuous-time simple random walks on a complete graph with 
pN vertices, where p € (0, 1) is fixed. He looks at the symmetrised distribution as in (1.1) and adds 
an exclusion constraint: there is no collision of any two particles during the time interval [0,(3]. The 
combinatorial structure of this model enabled him to express the free energy in terms of a cleverly 
chosen Markov process on No- Using Freidlin-Wentzell theory, he derives an explicit formula for 
the large-iV asymptotic of the free energy; in particular he obtains a phase-transition, called Bose- 
Einstein- condensation, for large (3 and sufficiently large p. Toth's work inspired partly our approach 
and future work will be contributed to question of application of our main results to this setting and 
the large deviations for appearance of long cycles. Large deviations for integer partitions and cycle 
structures, where the random walk bridges with different time horizon are weighted, are obtained in 
[A06]. 

The structure of the paper is as follows. In Section 2 we present all our results. In Subsection 2.1 
we describe our main large deviation results and a couple of conclusions. The application to general 
quantum spin models and to a non-commutative version of Varadhan's Lemma is given in Subsec- 
tion 2.2. In Subsection 2.3 we study the special case for simple random walks on a finite set and 
in Subsection 2.4 we provide some basic facts about the space Dp and large deviations theory. The 
Section 3 is devoted to the proofs of our main results. In the appendix in Section 4 we prove a lemma 
on pair probability measures and an entropy estimation which we use in our proofs. 

2. The Results 

In this section we are going to formulate our main results. In Subsection 2.1 we present the large 
deviations results and important conclusions. Following in Subsection 2.2 we apply this to a non- 
commutative version of Varadhan's Lemma and give an application of our large deviations result for 
quantum spin models and in Subsection 2.3 we study the very important case when for finite time the 
large time rate function, the Donsker- Varadhan rate function, is the rate function for the large iV-limit 
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under the symmetrised distribution. At the end in Subsection 2.4 we provide some preliminaries about 
the topology in Dp and some notion on large deviations theory. 

2.1 Large deviations for symmetrised distributions We fix throughout the paper (3 > 0. We 
equip the space Dp with the Skorokhod metric. Our main aim is to encode the combinatorics for the 
sum over permutations for the symmetrised measure (1.1) with a sum over pair probability measures 
with equal marginals. In order to formulate the large deviations results we introduce the following 
notations. By *p(Z d x Z d ) be denote the set of pair probability measures on Z d x Z d and we let 

<P(Z d x Z d ) = {Q G ^{Z d x Z d ) : Q w = Q (2) } 

be the set of pair probability measures on 7L d x 7h d with equal first and second marginal, respectively 
QW(x) = ^2 yeZ dQ(x,y),x G z<i > and Q (2) iv) = ^2 xe zdQ( x ,y),y G Zd - The relative entropy of the 
pair probability measure Q G *P(Z d x 7L d ) with respect to the product Q w <8) m is given by 



Note that Q *— ► H{Q\Q^®m) is strictly convex. All the rate functions of our large deviations principles 
include this relative entropy of pair probability measures as the part coming from the combinatorics 
of the symmetrised measure Pj^g'. The other part of the rate functions comes into play in each of 
the principles from large deviations principles for a product of not necessarily identical distributed 
objects. Thus the encoding of the sum over permutations with a sum over pair probability measures 
represents a certain two-level large deviations principle. This is seen in the definition (1.1) of the 
symmetrised measure P^/j j where permutations are sampled uniformly and for each permutation 
there is a product of not necessarily identical distributions of single random walks with initial and 
terminal condition. 

On the level of path measures we define the following functional on the space of probability measures 
on the set Dp of path as 

J< 8ym) (//) = >f {H(Q\QV®m) + I™(ri} for vey(Dp), (2.2) 

where the functional Ip Q) is given by 

4 0) (M)= ™p {(F,fi)- Q{x,y)]ogVg >y (e< F M)} for \i G tyiDp), (2.3) 

where we write £ for and (F, <5 € ) = f D F(u)8^(du) = F(£). Clearly, li Q) is a Legendre-Fenchel 
transform, but not one of a logarithmic moment generating function of a random variable, hence 
there seems to be no way to represent this functional as the relative entropy of fj, with respect to any 
measure. I^ Q> and J^ sym) are nonnegative, and I^ Q) is convex as a supremum of linear functions. 

Let tt s : Dp — ► M. d be the projection tt s (uj) = uj s for any s G [0,(3] and u> G Dp. We denote 
the marginal measure of /x G ^3(.D/j) on Z rf by fi s = jj, o 7T.T 1 G *p(Z rf ), and analogously we write 
/xo a = n o (tto,^) -1 G *p(Z d x Z rf ) for the joint distribution of the initial and the terminal point of 
a random process with distribution fi. If we restrict the supremum in (2.3) over all F G Cb(Dp) to 
all functions of the form u i— ► g(u)Q,ujp) with g G Cb(^ d x R d ), we see that Q = /io,/3 if Ip Q \^) < °°- 
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Indeed, 



sup 1 ^2 9{x,y){^{x,y)-Q{x,y))\ 

g£C h (R d xR d ) l 



= sup { 9{x,y)^{x,y)- £ Q(x, y) logE X)W (e»« '^) } 

which implies that //o,/3 = Q- Therefore the infimum in (2.2) is uniquely attained at this pair probability 
measure Q, i.e. 

H(no,p\no ®m) + _ sup_ (n,F- logE£ 0>7r/3 (e F ^)) if = M/?, 

+oo otherwise. 



jta™>( M ) = J F & c h (D,) x ' ' ' (2.4) 



In particular, J^ sym) is convex. 

Theorem 2.1 (LDP for the mean of path measures Ln). Under the symmetrised measure P^Tg 
the empirical path measures (£jv)iv>i satisfy a large deviations principle on ty(Dp) with speed N and 
rate function J^ sym) . 

The proof of Theorem 2.1 is in Section 3.1. The proof does not rely on any Markov property of the 
N random processes, hence this assumption can be dropped. 

We also have a large deviations principle for the mean path level. Note that any mean path is an 
element in the space Dp([O,0\;R d ) due to the averaging of paths with values in the lattice For 
the path level we consider the continuous embedding of the space L>([0, /?]; M rf ) into L 2 ([0, 0\; R d ) (see 
[Dor96, Lemma 2.3] for details). The scalar product for the latter space is defined by 

{£,u)= / d S (£(s), W (*)>K« ( 2 -5) 
J o 

for [0,(3] — > R d , where (■, -) K d is the scalar product on M. d . In the following we also write £ s for 
£(s). We define the following functional on L 2 ([0, /3];M. d ) as 

7L sym V)= inf {H{Q\QV®m)+l£\u)\ for uj G L 2 ([0, (3];M. d ), (2.6) 
Qe<p(z d xz d ) p > 



where 



If\u)= sup {(M- £ Q(x,y)]og^Je^)} for u G L 2 ([0, 0\- R d ). (2.7) 



Then the large deviations principle for the mean path reads as 

Theorem 2.2 (LDP for the mean of paths Y/v). Under the symmetrised measure P^Tg' the mean 

(Yn)n>i of the paths satisfies a large deviations principle on L 2 ([0, (3];M. d ) with speed N and rate 

r (sj 



function J (sy '" 



The proof of Theorem 2.2 is in Subsection 3.2. The contraction principle [DZ98, Th. 4.2.1] yields a 
large deviations principle for the mean path from the one for the mean of path measures. However, 
the identification of that rate function from the contraction principle with the one in Theorem 2.2 
seems to be rather difficult task from a technical point of view. Luckily, our proof of Theorem 2.1 is 
so general that it can be slightly modified to give the proof for the large deviations principle for the 
mean path. For details see Subsection 3.2. 
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Denote by S(Z rf ) all bounded function /: 7L d — > R. On the level of probability measures on TL d we 
define the functional J^ sym) on the set of probability measures on 7L d as 

J ( ^\p) = inf {#(Q|Q (1) ®m) + J< Q) (p)| forpG*p(Z d ), (2.8) 
Qe«p(Z d xZ d ) L J 

where 

j£ 0) (p) = supj^/(^)- ^ Q(x,y)logE^(eW^)}. (2.9) 

feB(Z d ) x£Zd x ,yeZ d 

Theorem 2.3 (LDP for the mean of normalised occupation local times Z^). Under the 
symmetrised measure Pj^' the mean {Zm)n>i of the normalised occupation measures satisfy a large 
deviations principle on *$(Z d ) with speed N and rate function J^ sym) . 

The proof of Theorem 2.3 is in Subsection 3.3. Here, the same remarks as for the proof of Theo- 
rem 2.2 concerning the contraction principle apply. For details see Subsection 3.3. 

In the following remark we compare the symmetrised distribution with the i.i.d. case. 

Remark 2.4. For i.i.d. random walks with initial distribution m, the empirical path measure (Ln)n>i 
satisfies a large deviations principle with speed N and rate function 

IpM= sup {(F,/i)-log m(x)EUe F ®)\. 



This is an application of Cramer's theorem [DZ98, Th. 6.1.3] for the mean of the independent identically 
distributed random walks with initial distribution m. Note that Ip Q) > Ip, m for the pair measure Q 
defined as Q(x,y) = m(x)S x (y) for x,y G Z d , since 

- Y, Q(*,y)logE£>^))>-log £ m(*)E^(E^(e^))) 

x,ye1 d x,y& d 

= -log m(x)K^ x (e F ^). 

x,yeZ d 

In particular J^ sym) > Ip. 

Our large deviations Theorems 2.1-2.3 may be extended by considering a finite and positive measure 
m not necessarily normalised to one. However, more interesting is the question if we replace the 
conditional probability measure ¥x, y by the measure ^x,y{-) = P!r(-H{£/3 = y}) for any x,y G Z d . This 
is included in the following proposition. 

Proposition 2.5. Let m be a positive finite measure on 7L d and let g: l^xl^ R + be a bounded and 
continuous strictly positive function. Replace Fx, y by g(x, y)^x, y in the definition of the symmetrised 
measure P^f in (1.1). Then 

(i) Theorem 2.1 remains true with the rate function replaced by 



(ii) Theorem 2.2 remains true with the rate function replaced by 

uj^ inf {H{Q\Q^ ®m) + W\uj)- V Q(x,y)log g(x,y)\ . 
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(iii) Theorem 2. 3 remains true with the rate function replaced by 

>f {H(Q\Q W ®m) + jf\p)- V Q(x,y)\ogg(x,y)\. 

Proof. We will prove (i). Define the function F g (u>) = \ogg(uo,ujp) for any path u> G ^([0, 0\; TL d \ 
Then with probability one with respect to ®f =1 P^a^m) 

N N 

Y[g(xi,x ff{{) ) = exp (X^ogstfo .^ )) = exp (iV^, F 9 >) . 

i=l i=l 
Clearly, cj i— > F 9 (uj) is continuous and bounded. Hence the large deviations principle follows from 
[dHOO, Th. 111.17] for any Q G ^(Z d x The rate function follows as p ^ ij ym) (m) - which 
together with (2.4) gives the proof. The proof of (ii) and (iii) follows analogously with (i) and the 
proofs of Theorem 2.2 and Theorem 2.3 respectively. Note for (ii) that 

w ~^(4^ )= . { J ? ym)(/i) - w^iog^)} 

_inf {^ sym) M- V Q(x,y)log 5 (x,y)} for to G L 2 ([0, /?]; Z d ). 



where $ : ^(-D^) — > -D([0, /?]; M d ) is the continuous mapping for the contraction principle, compare the 
proof of Theorem 2.2 in Section 3.2. Note that we used the fact that Q = /io,/3 if ^ (/-0 < 00 • ^ 

From the previous proposition we get. 



Proposition 2.6. Let m be a positive finite measure on Z d and let g: M. d x M. d — > R + be a bounded 
and continuous strictly positive function. 

(i) 

11 N N 

o-e6Ar xi& d ,\<i<N i=t i=1 (2 11) 

= - >f \H{Q\Q^®m)- V Q(x,y)log</(x,y)). 

(ii) T/ie unique minimiser of the rate function fi il sym) (/x) — ^ x y6Zd (i 0j p(x, y) log g(x, y) is given 

fry 

m°= E ^(^i^/r 1 . ( 2 - 12 ) 

where Q° G ^(^ x Z d ) is the unique minimiser on the right hand side of (2.11). Under the 
symmetrised measure the sequence (Ljv)jy e N converges in distribution to the measure fi° 

as N — > oo. 

Proof. (i) Proposition 2.5 gives that the left hand side of (2.11) equals — inf^gfp^^ 1^ (fi) — 
{lM),p, log 5)- If we use (2-4) and substitute Q = fiQ^ we get 

- inf {h(Q\QV ® m) - V Q(x,y)log<7(x,y) 

+ mf sup ((/i,F-logE^ o _(e F ^)U. 
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The latter infimum over \x is equal to zero. To see this, pick F = to get the lower bound. To get the 
corresponding upper bound take /x = J2 X y& d Qi x i vWx,y ° £ _1 and use Jensen's inequality to get 



We are going to prove that fi° is the unique minimiser of the rate function jjl ^ J^ sym) (^) — 
yeZ d Ho : /3(x,y)log g(x,y). This proves then both (ii) and (hi). For that, let [i G ty(Dp) be a 
zero of li sym) . As the relative entropy has compact level sets, there is a Q° G ?($(Z d x Z d ) that 

minimises the formula on the right hand side of (2.11). As \n) < oo, we have /io,/3 = Q° and 
hence 

= /fV)= ^p {i^F-logE? (e F ® 
FeC h {D ) L 

Clearly, F = is optimal and the Euler-Lagrange equations yield, for any g G Cb(Dp), 

which identihes fi as /j, . □ 

2.2 Non-commutative Varadhan's Lemma with Bose-Einstein Statistics In this section 
we use our large deviations principle for the mean paths under the symmetrised distribution IP^f 
to derive a non-commutative version of Varadhan's Lemma with Bose-Einstein statistics. Let h be 
a self-adjoint n x n matrix with h XtV < for all x ^ y, x, y G G = {1, . . . , n}. We define a Markov 
process on the finite index set G with transition probabilities 

-hy iX 5t , if y ^ x 



analogously 



Pt- 



w+*> -.KM (2 . 13) 

,{x, y) = P«(0 = y|f («) = i) = (e-<'-> S ) for x, 9 € G, 

where the matrix h is defined by h y>x = h V)X for y ^ x and /i XjX = — X^rr hz,x- F° r later convenience 
we let A: [— n,n] -tta continuous function such that X(x) = X x for x £ G, and we denote by ho a 
continuous function ho{x) = h X)X — h XtX for each igG, 

We let ./V Markov processes . . . , £ (JV) with transition probabilities (2.13) and time horizon [0, (5} 
be given. Let ¥ x ]y denote the conditional probability measure with density e^o h D(^(s))ds starting in 
x G G conditioned to terminate in y G G. We write E x ]y for the expectation with respect to the bridge 
probability measure F x ]y . 

We derive a large deviations principle for the mean path under the symmetrised distribution 
(compare (1.1)) 

^r = ^E £■■■ £ (g)™(^ p £w 

where the initial distribution m is defined by m(x) = — for x G G. 

We consider mean-field type interactions for the N Markov processes . . . , £ (iV) of the form 
iV ^/(f Eili^'')^ f° r some bounded continuous function /: R — > R. In the following we write 
jg(sym,h) £ Qr ^ e eX p ec t a tion with respect to the symmetrised distribution. 

Theorem 2.7. Fix (3 > anc? n G N. Ze£ h be a selfadjoint n x n-matrix with h x ,y < /or x ^ y for 
all x,g G G. 
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(a) The mean paths, Yn, under the symmetrised distribution F^^ satisfy, as N — ► oo, a large 
deviations principle on L 2 ([0,/3];M) with rate function I^ ym,h) defined by 

I^ h \co)= inf {h(Q\Q w ®m)+A Q ' h) (Lo)\, (2.14) 
Qeq3(GxG) L > 

foruj £ L 2 ([0,f3];R), where 

If> h \u)= sup UgM- E Q{^y)\og^Uo 9{s)as)As\\ 

(b) 

lim L logE% y ^ h) (e N fott J kZli& ) ) ds ) = S up { f f(cj(s))ds - I ( ^ h) (cj)\. (2.15) 

N^ooN \ J weL2([0,/3]iR) 1 JO J 

Proof. (a) is a direct application of our main Theorem 2.2, and (b) is an application of Varadhan's 
Lemma [DZ98, Th. 4.3.1]. □ 

We outline how this large deviations principle gives a non-commutative version of Varadhan's Lemma 
under Bose-Einstein statistics. Let p be the state on the algebra M of all complex n x n matrices 
given by p(A) = Tr (e'^A) for A 6 M. for the given matrix h, and let Tre - ^ = 1. We fix a self- 
adjoint element x £ M and some continuous bounded function /: R — > R. This self-adjoint element 
x describes a mean-field interaction expressed through the mean matrix 

1 N 

where Xj, z = 1, . . . , N, is a copy of the matrix x. Further, let /i (JV) = X^i=i where hi,i = 1, . . . , N, 
is a copy of the matrix h. Hence /i (JV) and x (JV) both act on the iV-th tensor product of the n-dimensional 
single variable space. The symbol Tr + denotes the trace restricted to the subspace of all symmetric 
N- variables with respect to any permutation of their single indices. The restriction to these symmetric 
variables is called Bose-Einstein statistics. 

We shall calculate the trace via the Feynman-Kac formula and our previous results in Theorem 2.7. 
For this we need the following path measure 

which is the probability for the Markov process to start in x £ G and to terminate in y £ G. Note 
that this measure can be normalised with the function gp(x,y) = (e-^o h D(^(s))ds-^^ = y}) to 
obtain the conditional probability measure Fx]y. Here we apply our Proposition 2.5 in combination 
with Theorem 2.7. Note that this results in substituting the symmetrised distribution F ( ^p' h) with 
the symmetrised measure 

creSjvXieA x N £A i=l 

The Feynman-Kac formula gives 

l im _L l og Tr + (e-^ (iV) -^ (iV) )))= lim 1 log p^ (e"tf'<* E£l #)ds ) 

, rP , (2-18) 

= sup {/ /(^( S ))d S -/^ ) H + (Q,(/ /3 )}, 

weL 2 ([0,/3];R) L JO } 

where we replaced m in Theorem 2.7 by the counting measure Cou, and where (Q,gp) denotes the 
expectation of g@ with respect to the pair probability measure Q. 
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The analysis of the variational formula of the right hand side of (2.18) gives the following theorem. 
Theorem 2.8 (Mean-field interaction with Bose-Einstein statistics). 

lim ^logTr + (e-^ (iV) -^ (JV) ))) = sup { f f{u{s))ds - I^> h \u) + {Q,gA 

= /?sup jnf {f(u)-I ( Q ' h) (u)-H(Q\QV® Con)}, 
«GK Qey(GxG) L J 

where 

If h \u) = sup \au- V Q(x,y)log^ h (e a Io^l{^ = y})\ forueR. (2.20) 
aeR 1 x ,yeG J 

Trace formulas like (2.18) and (2.19) go back to the work of Cegla, Lewis and Raggio [CLR88] in 
which the authors use a combination of large deviations theory and group representation to derive a 
variational formula for the free energy of mean-field quantum spin systems. Inspired by their work, 
Petz, Raggio and Verbeure [PRV89] derived a non-commutative version of Varadhan's theorem using 
C*-algebraic methods. We thus have in our Theorem 2.7 and Theorem 2.8 derived a non-commutative 
version of Varadhan's Lemma and hence a variational formula for the free energy of mean-field quantum 
spin systems under symmetrised distributions, i.e. a version with Bose-Einstein statistics. 

Proof of Theorem 2.8. The proof follows from Theorem 2.7 and (2.19) and the analysis of the 
variational problems, which is done in the following Lemma 2.9 and Proposition 2.10. □ 

Lemma 2.9. Let Q G ^P(G x G). If ui € L 2 ([0, /?]; IR) is a constant function then the supremum in 
sup UgM - Y, Q(x,y)logE^ h (e-fo9(s)mds 1{ ^ = y} ) j 

ff e^([O,0;R) x,yeG 

is attained at a constant function g. 

Proof. Fix any Q G <p(G x G). Clearly 

g^ Yl Q(x,y)\ogE^ h (efo9(sms)<is 1{ ^ = y} ^ 

x,y£G 

is as a convex combination of logarithmic moment generating functions convex and continuous. We 
introduce the Haar basis {/ij}i>o for L 2 ([0, /?]; IR) consisting of the functions hi defined by ho(s) = 1 
for s G [0, 0\ and if 2 m < i < 2 m+1 - 1, 

r 2 m / 2 , if (3{i2~ m - 1) < s < p((i + \)2~ m - 1); 
hi{s) = l -2 m l 2 , if /3((i + i)2- m - 1) < s < /3((i + l)2~ m - 1); . 
l_ , otherwise 

Now for every e > there is a m £ N and a function g in the space H. m spanned by the basis functions 
/jo./ii,...,^-! such that If' h \u) < (g,co) -^ yeG Q(x,y)logEi' h (eIo 9(s)^s)ds 1{ ^ = + £ If 
uj(s) = u for all s G [0,/3] we get the Euler-Lagrange equations for the finite dimensional variational 
problem as 
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We observe that 

is symmetric if the function g is constant. This is easily seen by 

v Eg(-e^>l{^ = y» < ^ Eg(-e^l{^ = y» 

E "( e<9 ' ?> ) " x,fe- = v}) 

and a complementary lower bound with mm xye \Q(x,y) > (indices x and y with Q(x,y) = do 
not contribute at all). Hence, the equations in (2.21) for 1 < i < 2 m — 1 are trivially solved for any 
constant function g. The constant is determined through the first equation of (2.21). □ 

Proposition 2.10. The supremum in the variational formula (2.15) with the mean-field energy f is 
attained at a constant function uj(s) = mGR for all s G [0,/3]. 

Proof. We split the proof into two steps. In the second step we will show that 

I { * ym ' h \uj) > J I^ h) (Q s )ds, (2.22) 

where Q s (t) = uj(s),t G [0,(3], is the constant function with value co(s). 
Step 1 From (2.22) we get 

- /%M S ))d S + /^V)> [\-f(u(s)) + I^ h \Q s ))ds 
Jo Jo (2.23) 

inf ■ 



>/3inf {-f{u)+ljr- h) @)}, 



and hence 



sup { /%( W ( S ))d S -/J ym ' h) H}</3sup{/( U )-/^- h) (n)} (2.24) 



W eL 2 ([0,/3];IE 

Step 2 We now prove (2.22). We only need to show that 

lf\u)>J ip\Q a )ds (2.25) 

for any Q G *P(r 2 ). We fix Q G <P(r 2 ). For any u G R we get from Lemma 2.9 that 

4 Q) (n)= sup {an- V Q(x, y) logE^fe^^ll^ = y}) ). (2.26) 

The function ^ x ygG Q(x, y) log Eg ,/l (e a ^ ^ s - )ds ll{^ / 3 = y}) is clearly convex and infinitely differen- 
tiable with increasing first derivative. From this we conclude that there exists, for any e > 0, a 
function g G L°°([0, (3];R) such that 

I^(Q t ) < g(t)u{t) - Y, Q(x,y)logE^ h (e^Io = y}) + e for a.e.t G [0,/3]. (2.27) 

x,yeG 

As y G L 2 ([0,/3]; [0,m]) we get the lower bound 

If\u>) > {g,u) - Y Q(x,y)logE^ h (efo9(ms) d s 1{ ^ = y} y (2>28) 

x,y€G 

Therefore we are finished with our proof if we show that 

Ef' fe (ei'<ff( s ^) ds i{^ = y}) < J K^ h fe^^^) ds t{^ = y}^jdt. (2.29) 
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By continuity, we may assume that g 6 TL m for some m G N. Then g can be written as g = 
HlZi(hk,9)H[P(k - l)2- m ,(3k2~ m ]}. Hence, we need to show that 

2 m 

E^fVo' 3 9(*)i W ds i{^ = y }) < 2- m ^E£' h (e< hfc ' fl >^ (s)ds l{£ j g = y}). (2.30) 

k=i 

But (2.30) is an application of the Holder inequality for random walk expectations. Hence, (2.29) 
follows. Now, (2.29) implies, together with (2.27) and (2.28), the inequality (2.25), and finishes the 
proof. □ 

Example 2.11 (Quantum-Spin-1/2 variables and Telegraph process). We shall apply our 
results for the mean field free energy in Theorem 2.19 to the quantum spin 1/2 model introduced in 
[Dor96]. In [Dor96] the non-commutative central limit theorem for the following model was studied. 
We show in this example the extension to the non-commutative Varadhan Lemma to obtain the 
mean mean-field free energy. In the setting of Theorem 2.7 and Theorem 2.19 we consider the set 
G = { — 1, +1} of possible spin values and the following process on G with transition probabilities 

««* +«)-»kw-*)-{ ,_?«: Hii, 

equivalently 

Pt{x,y) = ^(1 + xye~ f ) ,x,yeG. 

From Theorem 2.7 and Proposition 2.10 we get for this process and any continuous bounded function 
/: [-1,1] ->R the result 

lim 1 log ^)(e^Sf =1 = gup f /* /Ms )) da - J<«™>)} 

N^oo iV ^ \ J we L 2 ([0,/3];R) WO ' (2.31) 

= sup{/(n)-/^ m ' t) (n)}, 

where 

lf^\u) = inf {H(Q\Q^ Cou) + sup {ua- V Q(x, y) logEf (e a ^^ As t{^ = y}) } }. 

(2.32) 

To analyse (2.32) we have to calculate the expectations Ex(e a fo €( s ) ds t{£p = y}) for any x,y G G. 
This can be done by simple matrix calculations as follows. 

E%(e a foZ( s ) ds K{^ = y }^j = (y^^-^-^^x), (2.33) 

where 

*' = ( o -°i ) and ax = ( 1 I 

are the two-dimensional Pauli matrices. The eigenvalues of the matrix (aa z — ^(1 — a x )) are A±i 



| + a 2 . In principle, direct calculations of all the expectations in (2.33) lead to an evaluation 
of (2.32). The resulting mean mean-field free energy is therefore given by 

lim 1 log ^'(e* = /?sup{/(u) - ^(1 - v 7 !-^ 2 )}, (2.34) 

compare [Dor06]. As an alternative to the tedious calculations we immediately get a lower bound via 
the following considerations. Let a G R and let u a denote the eigenvector for the eigenvalue Ai of the 
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matrix (aa z — \ (\ — a x )). Then 



o) 



" ' «a(6 



defines a martingale for the telegraph process. We insert in the expectation on the right hand 
side of (2.32) and obtain 

sup{/?(cm-Ai)- V Q(x,y)logE%(Mj?>l{& = y})}. 

As M x (Mp a) t{^p = y}) is a probability measure for all x G G due to the martingale property we can 
combine the logarithm term with the one from the entropy to get 

I% rm - t \u) = inf {(3(au-X 1 )+ V Q(x,y) log ). 

Hence, as the relative entropy of probability measures is positive, we get as a lower the Legendre- 
Fenchel transform of Ai, which is given by (3^(1 — Vl — u 2 ). 



2.3 A special case for the rate function In this subsection we study motivated by the results for 
the quantum spin models in the previous section a very important special case for the rate function for 
the mean of the normalised occupation local times when our N random processes are simple random 
walks on Z d whose generator is given by the discrete Laplacian. We let A C Z d be a finite box and we 
put m equal the counting measure Coua on A. Instead of the symmetrised probability measure P^Tef 
in (1.1) we consider the symmetrised measure 

ceSjviEieA x n g\ i=i 

Proposition 2.5(iii) gives a large deviations principle for the mean of the normalised occupation local 
times under the measure Man > i- e - we nave 

^ilOg(^0 Z - 1 ) = -^ ) ^( !) ), 

where 

4 y A m) b)= > f {H (Q\Q W ® Coua) + ■/£»}, p£?(^), (2.36) 

Qe«p(Z d xZ d ) 

with 

J { p Q {{p) = sup {/? V f(x)p(x) - V Q(x, 2 /)logE :c (ei'o ,3 ^) ds l{^ = y})). (2.37) 

Our goal is to identify the rate function in much easier and more familiar terms. It turns out that 
this rate function is the Donsker- Varadhan rate function I\ , defined as 

/ A (p) = ( * E r-^=i " ' ^ SUPP(P) C A , (2.38) 

[ +oo , otherwise 

where \x — y\ = maxi<j<d |xj — yi\,x,y G is the lattice distance. This is the rate function for the 
large deviations principle for the normalised occupation local time lp as (3 — > oo. More precisely, denote 
by £[o,/3] the path of the random walk and define the sub probability measure Px^ = P x (-|supp (lp) C 
A) = IPx( - |C[o,/3] c A). The the normalised occupation local time satisfies as — > oo a large deviations 
principle on ^P(A) with speed (3 and rate function I\—Ca [DV75-83], where C\ is just the normalisation 
inf pg tp( Z d) \ x,yeA: (y/p(x) — y/p(y)) for the rate function. 

\x-y\=l 
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The surprising result is that the rate function Jp y ^ for finite (3 but for large N under the sym- 
metrised measure equals the Donsker-Varadhan rate function. 

Theorem 2.12. Let A C 1 d be a finite set. Then j££°(p) = PI\(p) for any p G with 
supp (p) = A. 

For the proof of the theorem we need the following lemma. 
Lemma 2.13. Fix (3 G (0, oo) and a finite set A C Z d . Then, for all p G ^(Z d ) having support in A, 

4T\p)= { H(Q\Q™® Coua) + Jf{ (p)}, (2.39) 

w/iere 

4 Q) A (p)= sup {/3 ^ Q(x,y)logE a; (e/«f/fe) ds ll{e [ o^ C A}1{& = y})} . 

/ eIRA xGA x,2/€A 

Proof. Note that H(Q\Q W <g> Coua) = oo if the support of Q is not contained in A x A. Hence, in 
(2.36) we need to take the infimum over pair probability measures Q only on the set ^3(A x A). From 
an inspection of the right hand side of (2.37) it follows that the function (vector) / in the supremum 
can be taken arbitrarily negative outside A to approximate the supremum. Hence, we may add in the 
expectation the indicator on the event that the random walk does not leave A by time (5. But then 
the values of / outside A do not contribute. This shows that we need to consider only functions / 
that are defined on A; in other words, (2.39) holds. □ 

Proof of Theorem 2.12. We proceed in two steps. First we show that Jp Y ^\p) > PI\(p) for 
any p G *P(Z d ) with support in A. In a second step we show the complementary inequality for any 
p G <p(Z d ) with supp(p) = A. We start from (2.39). 

Denote by Aa the restriction of the discrete Laplacian to the set A with zero boundary conditions, 
i.e. for any / G M A , 

A A /(x)= Yl [f(*)-f(y)l xGA,/GE zd ,supp(/)cA. (2.40) 

x,y£A : 
\x-y\ = l 

For / G / 2 (A) = R A let Uf be the unique positive eigenfunction for the operator Aa + / for the 
eigenvalue A(/) given by 

A A (/)= sup <(A A + /)<?,<?) = - inf {\ Y, (g(x)-g(y)) 2 -(f,g) 2 }. 

q&l 2 (A): 1 1 «| 1 2=1 9 6!2(Z<i): 12 ^— ' > 

9t( W ||g||2 S upp(g)CA,|| s || 2 =l x,yezd : 

\x—y\=l 

„ (2-41) 

To check the positivity and uniqueness of the eigenfunction consider the operator Aa + ci + /, which 
is non-negative and irreducible for some c > 0. The theorem of Perron- Frobenius (see [Se81]) gives 
that Aa + cl + / has a unique simple eigenvalue and that the corresponding eigenfunction is positive 
in A and unique up to constants. Clearly, Aa + cl + f and Aa + / have the same eigenfunction, and 
the corresponding eigenvalues differ by the constant c. Hence, A(/) is simple and the eigenfunction 
Uf is unique and positive. 

We now introduce a martingale which permits a transformation of the random walk in the expec- 
tation of the variational formula in (2.39). The expression 

M W := e /,f /(6)d Se -/3A A (/) 1{e[M c A}^||| (2.42) 
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defines a martingale {Mj )p>o under P x for any x G A with respect to the canonical filtration (see 
the Markov process variant of [RY99, Prop. VIII. 3.1] or [KalOl] for martingale theory for Markov 
processes). 

We insert now M^ f) on the right hand side of (2.39), obtain an extra (3\^{f) and use the marginal 
property of the pair probability measure Q, which gives J2 X y Q( x iV) 1°S u f f (x) = ®~ Thus we see that 

J^l(p)= sup {/?(</,?> -A A (/))- £ Q(x,y)logE x {M^t{^ = y})}, 



x,y£A 

where (/, p) = YlxeA f(x)p(x) is the scalar product in A. Here E^ denotes expectation for the simple 
random walk with generator A starting at x. Note that by the martingale property of the 
measure ~E x (Mp f) l{£p = •}) is a probability measure on A for any x G A. Substituting this in (2.39) 
and recalling the definition of the relative entropy of the pair measures, we obtain that 

J^(P)= >f Bup{/3«/,p>-A A (/))+ £ Q(x,y) log jML Tv}- ( 2 - 43 ) 

The double sum in (2.43) is, because E x {M ( f) i{^ = •}) is a probability measure, an entropy between 
probability measures and therefore nonnegative. 

Therefore we get that 

4 sym) (p) > sup Uf,p) - A A (/)}. (2.44) 
/eR A L J 

Note that the map / i— > Aa(/) is the Legendre-Fenchel transform of 7a 5 as is seen from the Rayleigh- 
Ritz principle in (2.41). According to the Duality Lemma [DZ98, Lemma 4.5.8], the r.h.s. of (2.44) is 
therefore equal to @Ia(p) since it is equal to the Legendre-Fenchel transform of Aa- Hence, we have 
shown that J A aym) (p) > (3Ia(p) for any p G *p(Z d ) with support in A. 

In our second step we prove the complementary inequality. For that we construct a pair probability 
measure to get an upper for (2.43). The resulting upper bound is given by the supremum over any 
function / G M. A . Thus we are left to solve this variational problem for the upper bound. Let's turn 
to the details. 

For p G *p(Z d ) with supp (p) = A define the function u G R A by u(x) = \Jp{x) and the function 
/* G R A by 

f*(x) = r -V z for x G A. 

Then u = uj* is the unique positive eigenfunction for the operator Aa+/* with eigenvalue Aa(/*) = 0. 
Given these objects we define the function Q* : A x A — > [0, 1] by 

Q*(x, y) = u^H^E* (e# ^^{^ C A}1{& = y}) , x, y G A. (2.45) 

This function is obviously symmetric and sums up to one. By the martingale property of {M^ ')/3>o 
the first marginal is identified as 

Q* m (x) = U (x)E x (ef»f*^ ds l{t: m C A}«(^)) = n(x) 2 , x G A. (2.46) 



Therefore Q* G ^P(A x A), which gives the following upper bound for (2.43), when we apply the 



marginal property of Q* and separate the eigenvalue in the definition of Mi f \ The marginal property 



gives 

N D* ( nr ii\ lncr 

Q*W(x) 



\ i u(x)u(y) 



x,y£A 
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Hence 

J^\p)< sup \(3(f,p) + Q*(x,j/)log ) f[. 



(2.47) 



We will show that the variational problem on the right hand side of (2.47) is solved for / = /*. Due 
to the strict concavity we need to get the solution for the Euler-Lagrange equation, which reads 

E,((/ ( f v(t s )ds)etf f^ ds l{t m C A}1{& = y}) 

fj(v, P ) = £ Q'(*.y) v / f , WM u x ( 2 - 48 ) 

x^A E a (e/o '^itf^ C A}1{& = y}) 

for all v G R A . The right hand side of (2.48) for / = /* reads 



u(x)u(y)E x (( f v(Z s )ds)efo f^t{£ m C A}1{& = y}) 

= Tu(x) 2 E x (Mp f v(£ a )ds)= / d g r»(«(6)), 

v Jo 7 JO 



(2.49) 



where E (/ * ) is the expectation with respect to the transform with the martingale (Mg )/?>o (compare 
remark following Prop. VIII. 3. 9 in [RY99]), starting at its invariant measure u 2 (x) = p(x),x G A. 
Therefore E (r) (v(£ s )) = {v,p), because the transformed random walk does not leave the set A and is 
stationary when started in its invariant measure. Thus f = f* solves the variational problem on the 
right hand side of (2.48) and gives finally after short calculations 

. Aa/u 
u 



J ( ^\p)<P(r,P) = -P(^,u) = -^u(x)A A u(x) = p(2d- £ «(*)«(*/)) 
= 4 E (VvJx) - Vp~¥)Y = PUP) 



xeA x,yeA: 

lx ~ yl=1 (2.50) 



2 

| a:-i/ 1=1 

□ 

We will give a heuristic interpretation in terms of the cycle structure. The measure l~i^ y p in (2.35) 
admits a representation which goes back to Feynman 1953 [Fe53]; in fact he considered Brownian 
motions instead of random walks on Every permutation a G &n can be written as a concatenation 
of cycles. Given a cycle (i, a(i), a 2 (i), . . . , a k ~ 1 (i)) with a k (i) = i and precisely k distinct indices, the 
contribution coming from this cycle is independent of all the other indices. Furthermore, by the fact 
that Hxi,x a ii) is the conditional distribution given that the random walk ends in x a ^, this contribution 
(also executing the k integrals over x a iu\ G A for / = k — 1, k — 2, . . . , 0) turns the corresponding k 
random walk bridges of length [3 into one random walk bridge of length k(3, starting and ending in 
the same point x\ G A and visiting A at the times (3, 2(3, ... ,(k — 1)(3. Hence, 

«4e ®(//«r M 

' o-eejv fceN JA 

where /fc(cr) denotes the number of cycles in a of length precisely equal to k, and fix',y is the random 
walk bridge measure fix,y = Px( - 1{^ = u}) j restricted to the event P|z=i{6/3 e A}. (See [Gin71, 
Lemma 2.1] for related combinatorial considerations.) If /jv(f) = 1 (i- e -, if a is a cycle), then we are 
considering just one random walk bridge £ of length N/3, with uniform initial measure on A, on the 
event D^ife/3 e ^}. Furthermore, Y/v is equal to the normalised occupation measure of this random 
walk. For such a a, the limit TV — ► oo turns into a limit for diverging time, and the corresponding 



18 



STEFAN ADAMS AND TONY DORLAS 



large-deviation principle of Donsker and Varadhan formally applies. In [A06] large deviation results 
are obtained for the cycle representation for integer partitions. 

It has been proved by Siito [Sii93, Sii02] that the non-vanishing of the probability for long cycles is 
equivalent with Bose-Einstein condensation for non-interacting Bosons. [DMP05] studied the density 
of long cycles for mean-field interactions and [BCMP05] study large deviations for different cycles 
statistics in the grandcanonical ensemble. 

If a permutation a does not contain a cycle of length ~ N a , for some a > 0, its contribution is 
quantified with a different rate, compare [BCMP05] and for the canonical ensemble [A06]. In this 
way, Theorem 2.12 says that the large- N behaviour of /U^yv' ° Y^ 1 is predominantly determined by 
all those permutations consisting of just one cycle of length N. 

2.4 Preliminaries 

In this subsection we provide some definitions and basic properties. A function £ : [0, 0\ — > M. d is an 
element of Dp if and only for any to G (0, (3) 

lim£(i) = £(i ) and lim£(i) exists, (2.51) 
tJ.to t]t 

and, moreover, £(0) = lim^o and £(/?) = linn,^ £(i). Every function in Dp has at most a countable 
number of points where it is not continuous ([Par67]). On Dp one can introduce a metric, the 
Skorokhod metric, defined by 

p(S,?)= inf { sup |f(A(t))-£'(t)|+ sup |A(t)-t||, (2.52) 

where H is the space of increasing functions A : [0, f3] — ► [0, 0\. The Skorokhod topology is the topology 
induced on Dp by the metric p and turns Dp into a Polish space([Par67],[DS01]). One has the following 
characterisation of compact subsets of Dp. 

Lemma 2.14 ([Par67]). In order that a subset K C Dp be compact under the Skorokhod topology it 
is necessary and sufficient that K be bounded and closed and satisfy the condition 

lim sup ^((5) = 0, (2.53) 
5J.0 ^ K 

where uj^(5) is defined by 

ut(8) = sup |i£(t') - m a m") - m\\ + ^ P m - m\ 

t-S<t'<t<t"<t+6 1 J 0<t<5 



+ sup |£(t)-£(/J)|. 

f3-5<t<f3 



(2.54) 



For the convenience of our reader, we repeat the notion of a large-deviations principle and of the 
most important facts that are used in the present paper. See [DZ98] and [DS01] for a comprehensive 
treatment of this theory. 

Let X denote a topological vector space. A lower semi-continuous function / : X — > [0, oo] is called a 
rate function if / is not identical oo and has compact level sets, i.e. if 7 _1 ([0, c]) = {x G X: I(x) < c] 
is compact for any c > 0. A sequence (A^jveN of ^-valued random variables Xn satisfies the 
large- deviation upper bound with speed ajv and rate function / if, for any closed subset F of X, 

limsup — logP(AAT € F) < - inf I(x), (2.55) 

and it satisfies the large- deviation lower bound if, for any open subset G of X, 

lim inf — log F(X N G G) < - inf I(x) . (2.56) 
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If both, upper and lower bound, are satisfied, one says that (Xn)n satisfies a large- deviation principle. 
The principle is called weak if the upper bound in (2.55) holds only for compact sets F. A weak 
principle can be strengthened to a full one by showing that the sequence of distributions of is 
exponentially tight, i.e. if for any L > there is a compact subset Kl of X such that F(Xn G KfJ < 
e~ LN for any TV G N. 

One of the most important conclusions from a large deviation principle is Varadhan's Lemma, which 
says that, for any bounded and continuous function F: X — ► R, 



One standard situation in which a large deviation principle holds is the case where P is a probability 
measure, and Xjy = + • • • + Yn) is the mean of N i.i.d. ^-valued random variables Y\ whose 

moment generating function 



is finite for all elements F of the topological dual space X* of X. In this case, the abstract Cramer the- 
orem provides a weak large deviation principle for (Xjy)Nm with rate function equal to the Legendre- 
Fenchel transform of logM, i.e. I(x) = sup FgA >» (F(x) — logM(F)). An extension to independent, 
but not necessarily identically distributed random variables is provided by the abstract Gartner-Ellis 
theorem. 

For our main theorems Theorem 2.1, Theorem 2.2 and Theorem 2.3 we shall rely on the following 
conventions. For the random probability measure in y$(Dp) we conceive ty(Dp) as a closed convex 
subset of M(Dp), the space of finite signed Borel measures on Dp. M(Dp) is a topological Hausdorff 
space, whose topology is induced by the set C\>(Dp) of continuous bounded functions on Dp. Here 
Cb(Dp) is the topological dual of Ai(Dp). The set ty(Dp) of probability measures on Dp inherits 
its topology from A4(Dp). When we speak of a large deviations principle of ^(-D^-valued random 
variables, then we mean a principle on Ai(Dp) with a rate function that is tacitly extended from ty(Dp) 
to M(Dp) with the value +oo. Thus, in the variational formula (2.3) the Legendre-Fenchel transform 
with respect to the set C\,(Dp) appeared. For the mean path we embed the space Dp([0, 0\; M d ) 
continuously into L 2 ([0, (3];R d ), hence the dual pairing is here given by the L 2 -scalar product. 



In this section we prove our main Theorems 2.1, 2.2 and 2.3. In the first subsection we prove the 
large deviations principle in Theorem 2.1 for the empirical path measures. In Subsection 3.2 we prove 
Theorem 2.2, in Subsection 3.3 we prove Theorem 2.3 both on the basis of the proof for Theorem 2.1 
and the contraction principle. The last subsection, Subsection 3.4, is devoted to the proof of the 
exponential tightness of the distribution of the empirical path measure under the symmetrised measure 
(1.1) and the exponential tightness for certain products of not necessarily identical distributed objects 
coming form the two- level large deviation method occurring in the proof of our main Theorems 2.1 
and 2.3. 

3.1 Proof of Theorem 2.1 

For the proof of Theorem 2.1 we have to show the following inequalities and properties, 
(1) For any open set A C ty(Dp), 





3. Proofs 



liminf log¥% y ^(L N € A) > - inf /< s: 




(3.1) 
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(2) for any compact set F C y$(Dp), 

lim sup log F% y ™ } (L N G F) < - inf /J ym) (p) . (3.2) 

(3) The sequence of probability measures (P^ r y ^' ) o Zy^)jvgf^ is exponentially tight. 

The symmetrised measure IP^T^ is not a product of i.i.d. measures, hence we first transform 
the probability measure into a probability measure which permits an application of a version of the 
Gartner-Ellis theorem. This will be done in Section 3.1.1 for the lower bound and in Section 3.1.2 
for the upper bound respectively. The technique applied here is a purely combinatoric one and it is 
inspired by the works of [Gin71] and [T6th90]. The counting arguments for the proposed combinatoric 
methods are given by the results in [A01], where combinatoric counting was used to evaluate partition 
functions in the microcanonical ensemble. For the whole proof of Theorem 2.1 the initial distribution 
measure m G ^3(Z d ) allows some kind of compactification. The exponential tightness will be proved 
in Subsection 3.4. 

3.1.1 Proof of the lower bound of Theorem 2.1. 

Fix an open set A G ty(Dp) and a finite set A C Z d with m(A) > 0. We get the lower bound 

1 N N 

F%^(L N eA)>m(A) N - £ £ •■■ E II^A(^)((8)Pl, CTW )(^GA), (3.3) 

ctGGjv ^iGA a;jvGAi=l i=l 

where m\(x) := m(x)/m(A) for any x G A defines a probability measure on A. The terminal points 
x a u\ are in A as well for any a G &n, %i £ A, 1 < i < N . We want to employ the following combinatoric 
scheme. We shall rewrite the sum over permutations with a sum on pair probability measures, where 
we are only asking for the frequency of transitions from an initial point in A to a terminal point in 
A. This frequency can be expressed in an easy way by a pair probability measure, i.e. a probability 
measure on A x A. Let 

= £(A x A) n 1< XA , (3.4) 

be the set of pair probability measures Q with equal marginals such that NQ(x,y) G No for any 
Q G ^P^' and any x,y G A. Note here, that at most entries of (Q(x, y)) XjJ/6 A are nonzero for any 
Q G y$p. Further, for a configuration x = (x\, . . . ,xn) G A n and given pair probability measure 

Qe&P, let 

& N (x, Q) = {a G & N - Vx, y G A : |J{z : x; = x, x a ^ = y} = NQ(x, y)}, (3.5) 

be the set of permutations which are admissible with a given configuration x and pair probability 
measure Q. Define the distribution 

»^-(®KP M ) <"> 

x,yeA 

for any pair probability measure Q G <P^. Note that 

N 

(g)P^ CT(i) = Pg iJV , Q G SjftV € ©at(x, Q),x e A N . (3.7) 
1=1 

Clearly the measure in (3.7) does not depend on a G ©at as long as a G &n(x, Q). Furthermore Q (1) 
is the empirical measure of the configuration x G A^, hence 

N 

il m ^) = U m ^ NQ(1){x) - 

i=l x£A 
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We insert a sum over Q G ^3^° in (3.3) and continue the estimation in (3.3) 
= m(A)" £ II-a^"^^^) E 

n m (JV) XGA SjGA, 
Vt+>A l<i<JV 

(3.8) 

The measure Pg ^ is a mixed product of single distributions of the random walks with given starting 

and terminal points. Recall that W% >y is the conditional probability measure for a random walk starting 
in x with terminal location y. Now we need to compute the counting term 

E $& N {x,Q) 
N\ 

l<i<JV 

because the remaining factors in (3.8) do not depend on x = (xi, . . . , xjv) G A^. We do this with the 
help of an additional sum over configurations uj G A N with x a = uj, where x a = (x<j(i), . . . , x^jv)) G 
(A) N for any permutation a G ©at. Thus 

E E 1= E E E % CT =^}(0-)%x,j/GA: #{*: a; i = a: , I1 ; i =y}=JVQ(x,j,)}(^) 

zGA^ ae6jv(S,Q) xeA^ wGA^ ctGBat 

= E EE % CT =^}( CJ ) ]1 {V 2; ,s/GA: x i =x,u H =y}=NQ(x,y)}{^) ^ ^ ^ 

~eA N : ueA N : <tG©jv 

L(x)=QW L(w)=QW 

where in the second line we took x G according to the marginal Q (1> , i.e. L(x) = Q w . Here, 
L(x) = Q w means that the first marginal of Q equals the empirical measure of the configuration x, 
i.e. 

1 N 

Q (1) (y) = ^E^(y) for y £ A - 

The same holds for the configurations uj G A^, i.e. the sum is restricted to L(uj) = Q w . Now we have 
to count the single terms in (3.9). We start with the sum over the permutations. For fixed x G A N and 
fixed uj G A N the number of permutations a G &n with x a = uj is exactly Y[ X £\(NQ^(x))\, which 
follows from the fact that for any x G A exactly (7VQ (1) (x)) times the x appears in the configuration 
x and uj giving for the first permutation NQ (1) (x) possibilities, for the second one iV(5 (1) (x) — 1 and 
so on. The sum over the configurations x G A N with Q (1> = L(x) is just the multinomial distribution 

N \ N\ 



K NQWj' U(NQW{x))V 

xeA 

where YlxeA NQ^(x) = N. In (3.9) the sum over the configurations x G A N gives the multinomial 
factor, because all the other terms do not depend on the configurations x anymore. Thus, together 
with the counting over the permutations we get 

E E 1 = N[ E 1 {^x,yeA:#{i:x i =x,a Ji =y}=NQ(x,y)}(^)- (3.10) 
*i£ A > a£& N (x,Q) ^eAN 

It remains to count the configurations uj relative to the given configuration x and frequency NQ(x, y) 
for any x, y G A. This number corresponds to the number of Euler trails in a complete graph or the 
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number of configurations in a microcanonical ensemble specified through the pair measure Q (cf. [A01] 
and references therein), and equals 

2^ H\/x,yeA: x i =x,u li =y}=NQ{x,y)}\^) ~ jj eA {NQ(x,y))\' ^ ' 

L(u) — x 

Thus we get from (3.8) 



F^\L N eA) 

>m(\\ N V" nm.M^'W UxeA( N Q W ( X )) [ p/3 ,r ^ (3.12) 

>m(A) ^ H m A(x) ^ = ^]vO?^V wG } ' 



Applying the Stirling formula for A! ~ N N+1 / 2 e N y/27r, we get after some computation and estimation 

n rn J T \NQW(x) Uxe\( N Q W ( x )) 1 > p C\A\ 2 logN -NH(Q\QW®m A ) 



x£A 

where 



H{Q\Q^®m h )= Y, Q(^y)log= Q(r " /; 



is the relative entropy of the pair probability measure Q G *p(Ax A) relative to the probability measure 
Q w (g) ?7iA G *P(A x A), and where C > is an absolute constant. Hence, 

^^L^ GA^m^e"^ 210 ^ J] e- NH WQ (1) ® m ^ Q N (L N e A). (3.13) 



In (3.13) there is a sum over all pair measures in ■ For any pair measure Q G ty{ r L d x Z d ) with 
finite relative entropy -ff (Q\Q®m) and finite functional we need an approximation with a sequence 
(Qn) JVeN of pair probability measures Qtv G ^Pa converging to Q as A — > oo in an appropriate way. 
Additionally, we let (ATv)ArgN be a sequence of sub boxes Aat C Z d such that An f Z d as A — ► oo. 

This will be done in the proof of the following proposition. 

Proposition 3.1. Let (An)n^n be a sequence of sub boxes An C Z d such that Aa? | Z d as N — ► oo 
and such that there is e > with |Aat| 3 /A = N~ e and 4|Aat| 2 /A < \ for all A G N. Then for any 
open set A C ^(Dp) we have 

liminf-J-log Y e- NH W QW ® mA N^ N {L N G A) > - inf inf {h(Q\Q w m) + 



Qe<p A 



(JV) 



N 



(3.14) 



Proof. Let A C y$(Dp) be an open set, and let minA and Q G ^P(Z d x Z d ) be given. In the case 

;/ii/V: ' '" ; 1 >P 



H(Q\Q (1) 0m) = +oo or Jg (/x) = +oo the assertion (3.14) follows immediately. We assume therefore 
that the relative entropy H(Q\Q W <g) m) and the functional Jg Q) (^) are finite. The first assumption 
implies that supp(Q) C (supp (m) x supp(m)), because 

H(Q\Q W ® m) = H(Q\Q W ® Q (1) ) + H(Q w \m), (3.15) 

where ff(Q (1) |m) = Yl x & d Q W ( X ) l°g • Hence, Q £ ^(supp(m) x supp(m)). 

Our strategy is as follows: 
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Step 1: In the first step we construct a sequence (Q N )NeN of pair probability measures Q N G , 
/at 

Step 2: lower bound 



such that Qffi — > Q weakly in the sense of probability measures as iV — > oo and such that the 



Ijminf llogP^,^^ e A) > (3-16) 

for any fi £ A follows from a lower bound of a large deviations principle for the distribution of Ln 
under the measure F 13 , m , and that 

Step 3: 

liminf -j- log f e -^ (Q ^ )|Q ^' (1)) ® mA iv )N ) > -H{Q\Q W ® m). (3.17) 

N— »<X> N \ ) 

We will start with step 1. 

Step 1: Construction of pair probability measures Q^ 1 ' G ^Pa'n 

This construction is splitted into two parts. First we construct for a given pair probability measure 
on Z ' x U 1 with equal marginals a pair probability measure on An x An with equal marginals such 
that it converges weakly to the given pair probability measure on Z d x Z d . This is done in Lemma 3.2. 
Ones we have a sequence of pair probability measures on A^v x Ajy with equal marginals we construct 
for this sequence a second sequence also with equal marginals but such that their single entries are 
integers when multiplied by N. 

Lemma 3.2 (Marginal construction). Let Q G ( ty{'L d x Z d ), A^v C Z d and xq G An and (tjn) a 
sequence in (0, 1) with t]n —> as N — > oo. Define the \An\ 2 entries of a function Qn- A-n x A 
[0,1] as 



N 



Qn(x,v) := < 



Then the following holds 



Q{x,y) , if x,y £ A N \ {x } 

QW(x) - T,zeA N \{x } Qn{x, z), if x G A N \ {x },y = x 

Q w {y) - Y,zeA N \{x }QN(z,y) , if ye A N \{x },x = x 

1 - E (*, v )eA%: Qn{x,v) , if (x,y) = (x ,x ). 

(x,y)jt(x ,x ) 



(3.18) 



(i) Qn € ^(A N x A N ) 

(ii) Let (An)ngn be a sequence of boxes Ajy C Z d with Ajy j Z d as N — ► oo. Then Qn — ► Q 
strongly as N — ► oo. 

(hi) Qat^o^o) > /or a// G N. 

Proof. (i) Clearly Qjv(z, y) G [0, 1] for all x, y G A w \ {x }, and Qat(x, x ) = J2zeA c N u{x } Q( x > z ) G 
[0, 1] and Q N (x ,x) = *EzeA%u{x } Q( z i x ) G t ' !] for a11 x e A N \ {x }. Also 

Qn(x ,x ) = 1-( Yl Q( x >v)+ E E {Q(x,y) + Q(y,x))^j g [0,1] , 

^>J/GAjv\{xo} xeAjv\{xo} 2/6A^U{xo} 

because the terms in the brackets are bounded by ^2 X eZ d Q(x,y). All |Ajy| 2 entries of the function 
Qn sum up to one, because 

E Qn(x,v)= E E ( 2 Q (1) W- E + 

x,y€A N x,yeA N \{x } xeA N \{x } yeA N \{x } 

+ <2aK x o,^o) = 1- 
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For any x G Ajv \ {^0} we get 

Q<$(x) =Q N (x,x )+ E Q(x,y)=Q w (x)- E Q N (x,y) + E Q(z,y) 

yeAjv\{a;o} yeA N \{x } yeA N \{x } 

= QV(x) =Q (2 \x) =Q N (x ,x)+ Yl Qn(v,x) 

y£A N \{x } 

= Qn 

and 

Q% ) (x ) = Q N (x ,x )+ E Q(y, x o)= E (Q w (y)~ E 

S/GAiv\{a;o} y£A N \{x } zeA N \{x Q } 

+ Qn(xq,x ) 

= i-( E Q(?>v)+ E ( 2 Q (1) W- E {Q(x,y) + Q(y,x)))) 

x,y£A N \{x () } a;GAjv\{a;o} yeA N \{x } 

+ E E 

J/GAjv\{a;()} z£A N \{xo} 

= 1- e q (1) w = i- E q ( n\x) = Qn\xo). 

x£A N \{x {) } x£A N \{x } 

Hence Qat G ^P(Ajv x Ajv). 

(ii) First, a direct computation gives 

Qiv(xo,xo) = l-2 E Q W ( Z )+ E 

zeA N \{x } y,zeA N \{x } 

= E E Q w ^ + E Q(*>v) 

zeA%U{x } y€A N \{x } y,z€A N \{x } 

= E Q w ^~ E Q(y> z ) 

yeA^ N u{x } 



E Q( x >y)> 

y,zeA c N U{x } 



this gives 



\Qn{xo,x ) - q(x ,x )\ = E Q( x >y)> ( 3 - 19 ) 

and hence 

E \Q(x,y)-Q N (x,y)\ < E Q0«M/) + E E f) + Q(v> x )) 

x,yeZ d (a;,j/)e(A Jv xA J v) c bsAjv^xo} J/GA^ 

<2 J] Q(x,y). 

(x,j/)g(A jv xA J v) c 

The assertion then follows since for any e > there is Nq G N such that 

E |Q(^y)-Qiv(x,y)| <£ for N > N. 

x,yeZ d 

was arbitrary. 

(hi) If Qn(xq, xq) = 52 y , ze A c N u{xo} y ) < ^ ^ or a ^ N G ^ we mmti Ply an tne entries Qn(x, y) 

for (x,y) G (Ajv x Ajy) \ {xo,xq} with the factor a/v = (1 — ??jv)/(l — Qn(xo, xq))- The resulting 
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entries are then denoted by Q N (x,y). Clearly, Q N (x ,x ) = 1 - J2(x,y)&A%\{xo,x } Qn(x,v) > Vn and 
a at — ► as N — ► oo. As 

\Q(x,y)-Q N (x,y)\ <2 ^ Q(x,y) + (1 - aw), 

x,yGZ d (x,y)e(A J vxA JV ) c 

all requirements in (i)-(ii) are satisfied by Qn- □ 

We now construct pair probability measures Qffi G ^Pa^- We apply the previous Lemma 3.2 for the 
choice i]n = (|Ajv| — l) 2 /N. Fix iV e N and denote the pair probability measure Qn in Lemma 3.2 
simply by Q. Hence, we have the property that here is xq G Ajy such that 

Q(x ,x )> (|AjV j v ~ 1)2 . (3.20) 

Now we define | Ajv | 2 — | Ajv| = : v components of a vector in [0,1]^ which satisfy the conditions 
(l)-(4) of (4.1) in Lemma 4.1. These coordinates define then according to Lemma 4.1 uniquely a pair 
probability measure in ^3^2 wnen the single components are elements of {0, 1/N, . . . , (N — 1)/A, 1}. 
Now (|Ajv| — l) 2 components are defined as 



Q w { ^ y):= lNQpy)i |oi . a]] rye Ay X{V{)} (;!21) 



where |_a;J is the largest integer smaller or equal to x G R. The remaining |Ajv| — 1 coordinates are 
given as 

Q W(i)) (x) := , for aJl x<E A* \ {*<)}, (3.22) 

where xo G Ajy as in Lemma 3.2 and in (3.20). We check that the conditions (l)-(4) in (4.1) of 
Lemma 4.1 are satisfied, where we order the coordinates in an obvious way. Since Qjv' (1)) (x) < Q (1) (x) 
for all x G A^ \ {xo} we have Yl x eA N \{x } Q^N ,W \ X ) < 1) i- e - (1) of (4-1) is satisfied. For any 
x,y G M+ one has [x + y\ > [x\ + [yj- Therefore, Q (1) (x) > T, y eA N \{x } and ^(z) = 

<2 (2) (x) > E^eAivVlxo} ( 3(j'' :r ) for a11 a; G Aat \ {x } imply that 

Q£ ,(1)) (s)> E <$°(*>v) and Q™W> E 

?/eA]v\{xo} ?/eA]v\{xo} 

for all x G A^v \ {xo}. Hence, also (2) and (3) of (4.1) are satisfied. For x,y G Ajy \ {xo} we have 
Qtf\x,y) > Qn(x,v) — = Q(x,y) — j?- Notice further that from the proof of Lemma 3.2 and our 
assumption (3.20) we have 

Q N (x ,x ) = l- E 2 Q (1) (z)+ E Q(x,y) = Q (1) (x )- Q( x i x °) 

x€A N \{x } x,yeA N \{x } xeA N \{x } ^ 

= Q(x ,x )> (|AjV j v - 1)2 . 
Hence, an estimation and application of (3.23) gives 

£ (2Q™>(*)- Yl Q { N } ^y))< E ( 2 Q (1) M- E <?(*,y)) 

xeAjv\{a;o} y€A N \{x } x£A N \{x } y£A N \{x } 

+ " A "' - 1)2 < 1, 

A 
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and therefore (4) of (4.1). We define the remaining 2|Aat| — 1 entries of the function see the 

proof of Lemma 4.1, as follows 

Q% ) (x,x ) = Q%- m \z)- Qn\^v) 

y&A N \{xo} 

Qn\x q ,x) = Q^ :{1) \x)- ^2 Q ( N } (y^ x ) for x G Aat \ {x }, 

s/eAjv\{zo} 

Q%\x ,x ) = l-( E 2Q™>(*)- E Qn\x,v))- 

x6Ajv\{aio} x,yeA N \{xo} 

Then the components in (3.21) and in (3.24) define uniquely a pair probability measure G ^3^° 
for JV6N. 

Weak convergence of (Qj^jveN 
For x,y £ Ajy \ {a?o} we have 

|Q^(^y)-Q(^2/)|<^ 

and from this we get 

Q (x, x ) - Q(x, x ) < | Q^' (1)) (*) - QV (x) | + £ | Q(x, y) - Q™ (x,y)\ 



y&A N \{x } 

J/GAjv\{x } 



(3.25) 



and analogously 



Q^ ) (x ,x)-Q(x ,x) < ^ forany xG Ajv\{x }. (3.26) 



Moreover, with (3.19) we get 

Ow'^cxo) - Q(x ,x ) < |Q^°(xo,^o) - (5at(xo,x )| + \Qn(,x ,x ) - Q(x ,x )| 

E o(«.v). <3 - 27) 

(x,y)G(A Jv xA JV ) c 

and therefore 

\Q ( N\^y)-Q^y)\<^^+ E SO^f)- ( 3 - 28 ) 

x,y€A N (x,y)e(A N xA N ) c 

and the assertion follows. 

Step 2: Lower bound via Gartner-Ellis Theorem We are going to use the Gartner-Ellis Theorem 
to deduce that 

1 s , i i s4 logpe «« (iA ' 6 - 4> -'" < '' ) - 

For doing this, we evaluate first the logarithmic moment generating function for a given F 6 C^(Dp), 
where we take the dual pairing ^(D^) with C^Dp) ([DS01, Lemma 3.2.3]). 

A(F):= lim ^logE* (e"< F ' L ">) = lim -J-logf TT l£ v (e F K>)"°" Wtt 

x,yeA N 

= lim £ Q^(x,y)logE^(e^)) 



x,yGAjv 

lim £ QW(x, y )logE^(e^))= £ Q(x,y) logE^(e^)), 



A^oo 

x,y£A N x,y£Z d 
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because the limit and the sum can be interchanged, and recall that Qn — ► Q (weakly) in sense of 
probability measures as N — > oo. Hence, A(F) exists, and 

A(F)= Y, Q(x,y)logElJe F ^) for all F€C h (D p ). 

x,y& d 

Also, it is easily seen that A is lower semi continuous and Gateax differentiable. Thus, [DZ98, 4.5.27] 
together with the exponential tightness for the sequence (P ^ (JV) o Ljj 1 )^^, which is derived in 

Lemma 3.5, gives 

liminfP^ (L N €A)>-I$» (jjl), 

N->oo Q N ,n 

where 

If\»)= sup {(F,»)-A(F)}, ntWDp). 

Step 3: Estimation for the relative entropy We are going to show (3.17), i.e. 

liminf-i?(Q^ ) |Q^' (1)) ®m Ajv ) > -H{Q\Q^®m). 

N— >oo 

We are going to estimate the relative entropy of the pair probability measure Qffi with respect to 
(5tv' (1)) <8> rriN- We intend to bound this relative entropy from above by the relative entropy of the 
restriction Q\a n of Q onto Ajy x Ajy with respect to the measure Q]^ ® m N plus some error terms. 
The latter relative entropy is known to converge to the relative entropy of Q with respect to Q ® m, 
see [DS01, Lemma 4.4.15]. We do this most conveniently with the help of the entropy 

H(Q ( n } ) = - E Q ( N\x,y)logQ^(x,y) 

x,y€A N 

of the probability measure and a well-know entropy estimation, which we cite and prove in 

Lemma 4.2 in the appendix. We write 

H(Q% ) \Q%'™®m N ) = -H(Q™)- £ Q$\x, y) log Q$ W) (x)m N (y) 

x,y£A N 

<-H(Q\ An ) + \H(Q\ An )-H(Q^)\- J2 Q ( N ) (x,y)logQ^ m \x)^m N (y) 

2 

F 

N~ 



....29) 

x,y£A N 



<-h(q\ An )- E Q^,y)iogQ^ m \x)m N ( y ) + ^^\og4N, 

x,yeA N 

where we applied Lemma 4.2 with our assumption that 4 ^ < i, the estimate 



x,y£A N 



N — 2> 



E \QW(*,v)-QkA*,v)\<^Tp-, (3-30) 



which follows from (3. 25), (3. 26) and (3.27), and a look at our previous assumption (3.20). Hence the 
entropy difference is bounded by 4| A^v| 2 log4N/N. Note 

- E Q%\*,v)b&Q l N W \*)™N{y) = - E Oz,y)iogQ^ (1) ^W(y). 

x,y&K N x, y eA N , 
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Further 



- E QN ) ( x ^y) l °sQ { N il) \ x ) m N(y) < - E QM^iogQ^VWO/) 



x,y^ N , x,y£A N , 

Q w (*)>-k Q {1) (x)>-k 



<- E QM*,y)iogQ (1 V(*W(y) + | E QI^Wiog v J (1) K ) 

<- E Q|A JV (^2/)logQ£( : r)m J v(y) + ^^, 

x,yeA N 

(3.31) 

because of 



Together with (3.29) and (3.31) we get finally the desired upper bound 



H(Q^\Q^ m) ®m N ) < -H(Q\ An ) - E Q\A N (x,y)logQ\ ( £ N W) (x)m N (y) + 0{N £ ) 

x,y€A N (3.32) 
= H(Q\ AN \Q\^ N ®m N ) + 0(N £ ). 

Clearly the right hand side of (3.32) converges to H(Q\Q W ® m) as N — > oo. Hence the assertion 
(3.17) follows. 

Finish of the proof of the proposition 

Now, (3.16) and (3.17) show that the left hand side of (3.14) is not smaller than 

-H{Q\Q^®m)-lf\n). 

Since this is the case for any \i G A and any pair measure Q € *P(Z rf x Z d ), the proof is finished. 

□ 

Finish of the proof of the lower bound of Theorem 2.1 

Now, Proposition 3.1 gives the desired lower bound for Theorem 2.1 when we apply it to the derived 
inequality in (3.13). Note that 

lim inf 4 log ( m(A N ) N e- c ^ 2 l ° sN ) = 

due to our assumption that \An\ 3 /N = N~ £ . Thus, 




□ 

3.1.2 Proof of the upper bound of theorem 2.1. For the upper bound we start with a finite box A C Z d . 
Later we will perform the limit A f Z d . The main task is to estimate the probability of the two events 
and to apply the combinatoric scheme introduced in the proof of the lower bound. Here, we are faced 
with the problem that a random walk may start in A and terminate in A or the complement A c , or 
start outside A and terminate outside or inside the box A. All events have to be estimated. To start, 
fix a closed set F C V( D /3), a box A C Z d and e > 0. 
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We split our sum over the iV initial points of the random walks into a sum over the initial points in 
the box A and over the complement A c for each single random walk. Thus we write 

E • ■ E = E E ■■■ E • (3-33) 

xiGZ d x N e1 d ae{l,c} N zieA"! x N eA a N 

The estimate the probability P^ /3 (Lat G F) from above, we write the probability as the sum of the 
following two events. The first event is the probability with additional indicator that the number of 
random walks not starting in A is more than eN, whereas the second event is the complement, i.e. 
the same probability with indicator that the number of random walks starting in A is greater or equal 
than (1 — s)N. Thus formally 

E - E = E E - E + E E - E • 

xiez d x N el d ae{i,c} N , x 1 eA a i x N eA a N a e{i,c} N , x!eA a i x N &A a N 

tt{i: a i = c}> £ JV a i = l}>(l-s)JV 

The probability for the first event can be bounded from above by the factor 2 N m(A c ) £N . For the 
probability of the second event we have to sum over all possible subsets I C {1, . . . , N} with |7| > 
(1 — e)N, hence 

N N 

' creSjv a€{l,c} N X!£A a i x N eA a N i=\ i=\ 

< 2 N m(Ay N (3 . 35) 

N N 

+ ^7 E E E E E n^^Oc^n 

cre6 N ic{i,-,n} ae{l,c} N : Xi£A a l X N £A a N i=l i=l 
\I\>l-eN I= {; ; a . = 1} 

We consider only those random walks in the product of random walk measures in (3.35), which are 
conditioned to start and terminate in the given box A. The contributions of the remaining random 
walks are estimated from above by one. For any given permutation a G &n and subset I C {1, . . . , N} 
with |/| > (1 — e)N we define the subset I c = {i G {1, . . . , N} : a(i) G /}. This means that we replace 
the product measure ®i^=i ^^,^(0 by the product measure (^ie^n/ ^i,x a u) ■ To perform this, we 
have to replace the empirical path measure Ljy by the empirical path measure Lj^i, given as 

Ll ' nT= \i^T\^ 6 ^- 

From |/| > (1 — e)N we get that |/o-n/| > (1 — 2e)N. We need some technical preliminaries to estimate 
the error of this replacement. First note that Dp equipped with the Skorokhod topology is Polish 
[DS01]. Next we need a metric for the probability measures on the Polish space Dp. Recall the Levy 
metric d on the Polish space ^p(Dp) [DS01], defined for any two probability measures fi, v G y$(Dp) as 

d(/x, v) = inf {/x(r) < z^(r 5 ) + 6 and u(T) < fi(T s ) +6 for all T = T C Dp} (3.36) 

where V s = {/i G ^(Dp): dist(/i, F) < 5} is the closed ^-neighbourhood of F. (By dist(/i, A) = 
vcd v( zA d(/U, v) we denote the distance to a set A C ty(Dp)). Note, NLn(A) is the number of random 
walk paths in any closed set A C Dp and therefore 

NL N (F) < \I a nI\L IanI + 2eN, 

which implies Ln(F) < Lj an j(F) + 2e. Thus, d(L>N, Lj an j) < 2e, and we get using |7| > (1 — e)N 

N 

((g>*t^))^ €F) < ( (g) K,^)(Ln G F*). (3.37) 
i=i ieiani 



30 



STEFAN ADAMS AND TONY DORLAS 



We insert this in (3.35) and execute all the N — \I\ summations over A aj with j £ I since they do 
not contribute anymore. All these contributions are estimated from above by one, and hence we are 
left with the |/| sums over those Xi with i £ I, i.e. Xi £ A for i G I. Hence, we write x = (xj)j e j for a 
configuration i£ A' and get from (3.35) 



%'(L N 6 F) < 2 N m(A')' N 



o-e6iv /c{i,...,Af> xgA 1 iei ieiani 

\I\>(l-e)N 

In the following we put \I a fl I\ = n and recall (1 — 2e)N < n < N. We observe that the product 
measure does not depend on the full information of the permutation a G &n- It depends only on the 
frequencies of indices i £ I such that x% = x and x a ^ = y for any x, y G A. These frequencies can be 
expressed with some pair probability measure 

Q G <£ ( A n) := <£(A x A) n 1n AxA ; 

n 

i.e. for example there are nQ(x, y) indices i G I a r\ I such that Xi = x and x CT (j) = y. Hence, we need 
to count those permutations that satisfy the constraint and further that 

Q m = L(x) = ^ K and Q (2) = L(a(x)) = ^ £ 

Note that here the marginals and Q (2) are not necessarily equal. But as n > (1 — 2e)iV we can 
estimate the difference of these marginals with the metric d on ^3 (A x A). The marginal differs at 
most by 2e, i.e. d(Q (1) ,Q (2) ) < 2e. Define the set 

y ( A e) = {Q£ x A): d(Q w ,Q (2> ) < 2e}. 

We rewrite the right hand side of (3.38) with a sum on pair probability measures in the set an d 
an additional sum over the subsets I C I with \I\ = n > (1 — 2e)N . 

^% y ^(L N G F) < 2 N m(A c ) £N 

+ m E E En ™(*r Q<1)w ( (2) &r Q( *' y) ) (L n s f 2s ) n -(x.) 

' o-eejv /c/c{i,...,jv}, Q e m(™' £ ) ieA x,yeA iei\i 

n=\I\>(l-2e)N A 



< 1{I = J g ft J} 



zeA 7 o-e©jv 

< 2^m(A c ) £jv + £ £ ™(^) n ° (1){s) ( ® (^,,) nQ(X * y) ) {Ln G F 2e ) 

(l-2e)Af<n<A r Q g( p^).£ xGA z.yeA 

x E E JJ] & n(x,T,Q) JJ m(xi), 

7c/c{i, r .,AT}, xeA 1 ' j € A/ 

„=i/i 



(3.39) 



where we introduced the set 

&n{x,I,Q) = {a G 6at: Vx,y G A: G /: Xi = x,x a ^ = y} = nQ(x,y)} 

of permutations admissible with Q and the configuration x G A^ on the index set I C I. 

Counting: We estimate now the cardinality of the set 6jv(x, I, Q). We fix a configuration x £ A 1 
and an index set I C / C {1, . . . , N} with \I\ = n, (1 — 2e)N < n < N. We evaluate the cardinality of 
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the set &n(x,I,Q) with an additional sum over those configurations uj £ for which Lj{u) = Q (2 \ 
here Lj(uj) is the empirical measure of the configuration u. This gives 

$& N {x,I,Q)= ^ \v x , y eK: )){j e 7: Xi=x,Wi=y}=nQ(x,y)}^ 

^(-)=<3 (2) (3.40) 

X E Vi^w-Viei}^)- 
o-eSjv 

The last term is easily estimated as 

E Vi^w.Viei}^ = E |{<7:/-»/bijective:w i = a; ff(i) Vi€7}| 
ffee " /c/,|/|=|/| 

x |{<7: {l,...,JV}\7-.{l,...,JV}\/bijective} 



< 



7 ') IJ( n ^ (2) ( x )) ! (^- n ) ! - 



xgA 



The first term in (3.40) is given as in the lower bound via a counting of Euler trails ([A01]). Therefore 

v ii ~ (u) < n, eA w (i) (x))! 

n {Vx,yeA: ${iel: x i =x,uj i =y}=nQ{x,y)} V" > - ^ gA ( n Q( X) ' 

L(i„) = q( 2 ) 

Note that the two previous estimations do not depend on the configurations x € A 1 as long as 

Q(i) = L 7 ( 5 ) = I^ (J 

The number of all these configurations is clearly equal n!/ n^eA^^ 1 ^))- ^ e s Pbt the sum over the 
configurations x 6 A^ in the last line of (3.39) into a sum on (xj) ig j £ A J and ( x i) ie j\j £ A^, and we 
estimate the term Y\ ieI \jrn(xi) from above by one. Hence, the last line of (3.39) can be estimated as 



E E ^&n(x,I,Q) JI m(^)< tt{(xi) ie y€A 7 :L((x i ) i6 7) = (1) } 

:{i,^..,jv}, zeA / ' igjAj /c/c{i, L ..,JV}, 

,=|/| n=|/| 

'l\ EUa (x))! EUa (nQ (2) (x))! (N - n)\ 



n 



l=n 

< N 



Ylx,yeK(nQ{x,y))\ N\ ^ 

V 'N\ fl\ 2 n\(N-n)\ EUaW 2) (x))\ 



n\ 2 n, eA w 2) w)! < Y iAP + c7 c2g£j y n, eA Q (2) (x)" Q(2)w 



«y Ux,yeA(nQ(x,y)y. - \\ x , y eM x ^) nQ{x ' yy 

where we used that there is a constant C > such that (^) < e CeN for all n with (1 — 2e)N < n < N, 
and where we used Stirling's formula and some constant C. Inserting all this in the estimation (3.39) 
gives 

N 

P%%\L N eF)< N \M*+c e 2CeN Y e~ nH ^ m ® Q(2) % <n (L n e F 2e ) + 2 N m(A c y N . 

n>(l-2e)N Q£ yM 

(3.42) 

Large deviations: 
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We now show for A C Z d fixed and any e > that 

hmilogf Y e -nH(Q\m®Q^) p P ^ F 2s )\ 

Q ^ (3.43) 
<- inf inf \H{Q\m®Q < - 2) )+ sup \(F, + V Q(x, y) log Eg v fe F ^U ), 

where <#f = {Q G <p(Ax A) : d(Q«, Q< 2 >) < 2e} is the set of pair probability measure whose marginals 
differs by at most 2e, and where the probability measure Q G is trivially extended to Z rf x Z d 
by zero outside of A x A. To see this, consider the logarithmic moment generating function of the 
distribution of L n under the probability measure Pq n , 

AW)(*) = logEg n (e^ L «>) =ny j Q(x,y)logEi y (e F ^) (3.44) 

for any Q G ^P^ 1 ' £) and any F G C\,{Dp). Now let Q n G ^P^ 1,£) be maximal for the mapping Q ^ 
e -nH(Q\ m ®QW) F ^^ Ln G F 2 £ )_ Then ^ gince get $p(e) ig compact) there is a pair measure Q G $J}W 

with limn^oo Q n = Q weakly. Clearly the limit 

A«(F) = lim -AW)(F)= J] Q(x,y)logE£>^)) 

exists, and is lower semi continuous and Gateax differentiate. Now the Gartner-Ellis theorem yields 
that 

limsup-logP^jL^GF 2 ^- inf sup {<F,/i>+ ^ Q(x,y)\ogV% >y (e F ®)}, 
™ n " ^ 2e Fec h (D ) «■ ' J 

because we may assume that F 2e is compact. We can do so, because Lemma 3.4 shows that there is 
a sequence of compact sets Ml C ty(Dp) such that 



lim lim sup - log ( sup Pg n (L n G M£)) 



-oo. 



r i i ( 3 - 45 ) 

inf inf {H(Q\m®Q {2) ) + I ( 3 Q \fi)\\, 



The cardinality of the set *P^' £) is clearly polynomial in n, and by continuity of Q H(Q\m Q (2) ), 
the assertion of (3.43) follows. 

We are now in the position to perform the N — > oo limit for the upper bound in (3.42). We get 
limsup — \og¥^\L N G F) < - min { - 2Ce - log2 - elogm(A c ), 

inf inf 

" eF2£ Qe<p A 

where we recall 

4 0) (M)= sup J] Q(^y)logE^(e^))}. 

The proof of the upper bound (3.1.2) of Theorem 2.1 is finished when we replace A C 7L d and e by 
sequences (£7v)jv e N and (Ajv)jveN with £jv — ► and A^r \ U 1 as N — > oo such that |Aat| 2 /A = A _<5 
for some <5 > and such that eArlogm(A^) — > — oo as N — > oo and use the following lemma. Recall 
that mjv G ^P(Aat) is the restriction of the initial distribution on the set An- 

Lemma 3.3. Fix a closed set F C ty(Dp). Then for any sequence (e n );veN satisfying ejy — ► as 
N — ► oo and any sequence (Ajv)jv e N wzi/i A^r T ^ d as A — > oo, 
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liminf inf inf ( H(Q\m N ® Q (2y ) + A Q) (fi)}} 

JV— »oo fi£F 2e Qgm( £ ) I ^ J J 

(3.46) 

> inf _inf {h{Q\Q {1) ®m)+lf\Li)}. 
^ eF QG«P(Z d xZ d ) L > 

Proof. Clearly, tun — > m weakly as A ^ oo. Now we pick approximating sequences of Q's and /x's 
and employ compactness arguments. Thus, for any A £ N pick /xjv G F 2ejv and Qn £ ^Pa^ sucn that 
the sequences {H{Q^\m,N <8>Qjy ) + ^ Qjv> (^Af))AfGN converges to the right hand side of (3.46) and may 
therefore be assumed to be bounded. The sequence {H(Q^\mN))Nen is bounded because of 

H(Q N \m N Q$ = fl-(Qj^lmJv) + H(Q N \Q% ® Q$). 

As Qjy has support in we have H(Q^\mN) = H(Q < ^\m), and thus the sequence (Qj^)jveN is tight 
due to the fact that the level sets of the relative entropy are compact (see Lemma [DZ98, 6.2.12]). 
As d(Q$,Q$) < 2ejv -► as TV -► oo, also the sequence (Qjy)jveN is tight. By boundedness of 
the sequence (H(Qn\Q n ) <8> Qjv))weNj also the set P := {Qn- A G N} is tight. Hence, according to 
Prohorov's theorem we may assume that Qn — ► Q as TV — ► oo for some Q G <p(Z d x Z d ). Since both 
Q# and Q$ -► Q< 2 > as A -> oo, we get that Q G <p(Z d x Z d ). 

For C > sufficiently large, the sequence (nn)Nefi is contained in the set 

ipeWDp): mf N /^(^)<C}. 

Now it turns out that this set is relatively compact. We are going to prove this fact now. Note that 
it suffices to find a family of compact sets Ml C Dp, L > 0, such that 

lim inf inf F Q) = oo. 

L->ooQePM£ 

We prove this in the usual way with the exponential tightness and a lower bound for a large devia- 
tions principle. The sequence (Ln)n& is exponentially tight under the probability measure Pq n , 
uniformly in Q G P (see Lemma 3.5). Moreover, it is easy to see that it satisfies a large deviations 
principle with rate function 1^ '. Indeed, note that the logarithmic moment generating function of 

Ln under the probability measure Pq^ n is easily shown to converge towards the function 



Yl Q(x,y)log^ y (e F ^) ,F£C h (Dp), 



whose Legendre-Fenchel transform is The Gartner-Ellis theorem then provides the proof for the 
large deviations principle. For L G N, pick a compact set Ml C ty(JDp) such that 

F QnM Ln g M i) ^ e ~ NL for aQ y L,N eN, Qn G P. 
Now the lower bound in the mentioned large deviations principle gives us that 

inf inf JjL Q) > - liminf — log P^ N (L N G Mf) > L, 

implying that the sequence (hn)ngn is tight. Therefore, we may assume that =>- ji as A — > oo 
with some \i G F 1 . Since fiN G F 2eN for any A G N and since £tv — ► 0, we even have /i G F, because F 
is closed. To finish now the proof of the lemma we employ the representation of the relative entropy 
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(3.47) 



as a Legendre transform (see [DS01, Lemma 3.2.13]. This gives 
inf inf {H{Q\m N ®QW) + lf\n)\ 

>(#,Qn) -\og(&, m N ®Q™) + l<S>,ii N )- QN(x,y)\ogE^ y (e F ^, 

x,yeA N 

where g G Cb(A7v X A^r) and F G C\>{Dp) are arbitrary. Hence, we get 
liminf inf inf {H(Q\m N ® Q (2) ) + li Q) (fi)\ 

> ( 5 ,Q)-log(e^,Q«®m) + (F, M )- £ Q(x,y) logE^(e F ^). 

rr,j/GZ d 

Since this holds for any g G Ct>(AAr x Ajv) and any F G C\,(Dp), the left hand side of (3.47) is not 
smaller than H{Q\Q W <g> m) + lf\p). Therefore, 

l.h.s. of (3.47) > inf inf \ H(Q\Q ® m) + I ( g Q> (n)\ , (3.48) 

and the assertion of the lemma follows. □ 

3.2 Proof of Theorem 2.2 We prove Theorem 2.2 in the following. Denote by * the continuous 
mapping 

vf : «p(D^) - £>([0,/?];K d ),/z ^ *(//) = / w/x(do;). 

Note that Y N = W(L N ) and recall that Dp = Dp([0, /?]; Z d ). The contraction principle [DS01, 
Lemma 2.1.4] ensures a large deviations principle for Yjy under the symmetrised measure IP^T^' with 
the rate function 

I% ym \u) = _inf \H(Q\Q^ ®m) + K Q) (oo)}, (3.49) 
Qe<p(z d xz d ) L M J 

where 

^ Q V) = inf I^(fi). (3.50) 

*(/a)=w 

Therefore we need to show that J^ sym) = lj* ym \ and for that it suffices to show that I^ Q) = I^ Q) . If we 

consider the class of functions F G C^Dp) of the form Ff(u) = ds(uj s , f s ) R d for / G L 2 ([0, f3];R d ) 
we get for /i G ^P(-Dg) 

4 Q) (/x)> sup {/ M (du;) [ ds(u s ,f s ) Rd - V Q(x,y)logEj J/ (e^ ds ^^V)) 

= 5 Q) (^)). 

If we now take the infimum over all probability measures fj, G ty(Dp) with *(/x) = w we get that 
i^ Q) > J^ Q) . To prove the complementary bound i^ Q) < I^ Q) seems to cause major technical difficulties. 
We therefore proceed in an indirect way. We show that both i^ Q) and il Q) are the rate function for the 

same large deviations principle. In the proof of Proposition (3.1) we have shown that the empirical 

HQ 



path measures Ljy satisfies a large deviations principle with rate function 1^' under the measure 

^Qn n ' w here Qn G ^P^2 * s ^ e se Q uence °f P a i r probability measure from step 1 in the proof of 
Proposition 3.1. According to the contraction principle the sequence (Ijv)jveN satisfies, under the 
measure ^q n jv > a large deviations principle with rate function . 
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Now we show directly that the sequence (Yiv)iVeN under the measure ~Pq n ^ satisfies a large devia- 
tions principle with rate function Ip Q \ which finishes the proof. Recall that Qn G with Qn — ► Q 
as N — > oo. We perform this in the usual setting of the Gartner-Ellis theorem. Let / G L 2 ([0, /?]; 
the logarithmic moment generating function is then 



A J v(/) = logE^ Ar (e^))=log( [J ^(e^^)' 

^^j/gAjv 

= iV J] Q JV (x,y)logEf )J/ (ei'o ,3ds ^^V), 



and hence 



A(/) = ^lim o lA JV (/)= ^ Q(x,y) log Ej^e^^^V^. 

It is easily seen that A is lower semi continuous and Gateax differentiable. The Legendre-Fenchel 

7(0 



transform is equal to /i* 3 ' '. According to Lemma 3.5, the sequence (Y/v)atgn is exponentially tight 



under (^q n n)n£N- Hence, the Gartner-Ellis theorem finishes the proof. □ 

3.3 Proof of Theorem 2.3 In this subsection we prove Theorem 2.3. We denote by ir s : Dp — > 
7L d ,s G [0,(3], the canonical projection it s (lu) = uj s for lo G Then the mapping 

T: %(Dp) -+ m d ),V » Tfjj) = - j\sii o ttJ 1 

P Jo 

is continuous and Zjy = T(Ljy). Then a large deviations principle for the mean Zjy of occupation 
local times with rate function 

J^ m \p) = inf {H(Q\Q {1) ®m) + Jf\p)\, (3.51) 

where 

^» = ^ (3-52) 

is given via the contraction principle [DS01, Lemma 2.1.4]. Therefore we need to show that J^ sym) = 
Jp yTa) , and for that it suffices to show that J^ Q) = Jp Q) . We relax the set of functions over which we 



perform the supremum. For a fixed probability measure /x G ty{Dp) we consider bounded continuous 



functions F G C^(Dp) of the form F(u>) = 4 Jq dsf(uj s ) for any bounded function / G B(Z d ) and 



u G Dp. Then for Q G <p(Z d x Z d ) 

/<?>(//)> sup {/ /.(do;)! f P dsf(u s )- Q(x,y)logEUe? J ° m)ds )} 
= J«>(T( M )). 

If we now take the infimum over all probability measures G ^(D^) with T(/i) = p we get that 
J^ Q) > J^ Q) . To prove the complementary bound J^ q) < J^ Q) seems to cause major technical difficulties. 

We therefore proceed in an indirect way. We show that both J^ Q) and ji Q) are the rate function for the 
same large deviations principle. In the proof of Proposition (3.1) we have shown that the empirical 

HQ 



path measures satisfies a large deviations principle with rate function il" 5 ' under the measure 



Fq n N , where Qn G is the sequence of pair probability measure from step 1 in the proof of 

Proposition 3.1. According to the contraction principle the sequence (Zjy)iy € fq satisfies, under the 
measure ¥q n n , a large deviations principle with rate function J^ Q) . 
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Now we show directly that the sequence (Zn)n^ under the measure ^q n n satisfies a large devia- 
tions principle with rate function Jg , which finishes the proof. Recall that Qn G ^Pa^ with Qn ~~ * Q 
as N — > oo. We perform this in the usual setting of the Gartner-Ellis theorem. Let / G B(Z d ) any 
bounded function, the logarithmic moment generating function is then 



A J v(/) = log^ jViAr (e^))=log( [J Ky{^ m)dS ) ) 

x,y€A N 

= N ^ Qiv(x,y)logEj J/ (ei'o /3 ^) ds ), 



and hence 

A(/) = iv lim o lA JV (/)= Q(x,y)logE^eIo^ ds ). 
It is easily seen that A is lower semi continuous and Gateax differentiable. The Legendre-Fenchel 

t(Q 



transform is equal to j£ Q) . According to Lemma 3.5, the sequence (Zn)^^ is exponentially tight 



under (Fq^ ^Nm- Hence, the Gartner-Ellis theorem finishes the proof. □ 

3.4 Exponential tightness We prove in this subsection the exponential tightness of the distributions 
of the empirical path measure L jy under Pjy^ and under Pq n for any pair measure Q G • For 
the first result we use the compactification given by the initial distribution m of the random walks. 

Lemma 3.4. The empirical path measures (L^)]^^ are exponentially tight under the symmetrised 
measure P^p^ ■ 

Proof. For / G N, choose a subset A; C Z d such that m(A^) < e _i . Furthermore, choose Si > so 
small that 



sup Pf )W ( sup |Ct' — £t| A \£t" — £t| + sup |£(i) — £(0)| + sup m)~m\ > y) <e"' 2 
i.yGA; H~8 l <t'<t"<t+8 l 0<t<S t P-&i<t<P u 

where we write £ = for a single random walk. Consider the set 



1 



Ai = [u G Dp : uj(0) G A;, G A/, sup \u t > - u t \ A \u t n - u t \ < —, 

t-Si<t'<t"<t+S t <j£ 

sup |^)-w(0)|<^ sup \u(t)-Lu((3)\<±-}. 

0<t<5 t M l3-5i<t<f3 M 

According to a well-known characterisation of compact subsets in Dp, compare e.g. [Dor96, Lemma 2.1] 
or [Par67], A\ is relative compact in Dp with respect to the Skorokhod topology. Now put Ki := {// G 
?p(Dp) : < j} and note that Ki is closed by Portmanteau's theorem. Let M G N be given and 



consider Km '■= f] Ki- The set Km is tight and by Prohorov's theorem Km is compact. We shall 

l=M 

show that W%*p\L N G K C M ) < e~ MN . Observe that 



{L N G Kf} C | U{i G {l,...,iV}: £ w G A?} > 



AT' 

7. 

C {«{,: e« G A?} > ^} U {«{•: € A?} > f } U tf> G A,,$> G A,, 

sup |fr - 61 A |6» - 61 > ttj, sup |£(t) - £(0)| > -, 
t-Si<t'<t"<t+Si 61 0<t<8 t <J £ 



sup \t( t )-t(p)\>h>%}- 



p-8,<t<i3 
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Clearly, 



Pj^GKt : $> € A^} > f ) = ^(Hi: tf> £A?}>^ 



3r ^ vnL " u 1J " 31- 

N N 



^ E ^iE E ii^l^tVie^^Af) 



AT! 

/C{1,...,JV}: creSjv xe(z d ) N 1=1 1=1 



(3.53) 



(3.54) 



< J] m(Af)l 7 l < e~ lN / 3 2 N . 

Furthermore, 

^N?Mi-- $ e A,,^ e A/, sup |& - 61 a |6» - 61 > 4 

sup |£(i)-£(0)|>l sup |£(i)-£(/J)|>I 
0<t<<5 ; ol p-8 t <t<l3 3L 

1 W ^ 

|/|>|f ' aee N xe(z d ) N »=i »=i 
sup 16' -61 A |6" -61 > ^7, 

t-5|<f <t"<t+<5; <3' 

sup iew - ewi > tJt, sup iec<) - oti > t^t) 

0<t<5i 6L /3-Si<t</3 6U 

< Yl su p E II m ^ H F -^( su p i& - &i A i&" - 6i > tJt, 
sup \m - m\ > L sup m - m\ > ^) 

< e- lN / 3 < e- lN ' z 2 N . 

Hence, 

oo oo 

V%p\L N € K C M ) < F N,p\ L N e Kf) < 32 N J2 e" W/3 < 6 x 2 N e~ NM / 3 < e~ NM / 5 

l=M l=M 

for all large N if M > 24. This ends the proof. □ 

The following exponential tightness is due to the product structure of the probability measure Pg N 
for any Q € • Here we have a product of not necessarily identical distributed random walks. 

Lemma 3.5. Let (Qjv)jveN be any sequence of pair measures Qn & ^Pa^ on An x An with An T ^ d 
as N — ► oo. TTien i/ie family of distributions of Ln (respectively of Yn and of Zn) under Pq n N is 
exponentially tight. 



Proof. We will prove the case for the empirical path measure Ln- The proofs for Yn and Zn follow 
analogously. Our proof is an adaptation of the proof for the i.i.d. case (see [DZ98, Lemma 6.2.6]). 
From the previous Lemma 3.4 we have a compact set A\ C Dp and a subset A^ C Z d such that 

sup P^(£ G ^) < e" 2Z V - 1). (3.55) 
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The set M t = {v G *#(Dp): v(Af) < 1/1} is closed by Portmanteau's theorem. For M £ N define 
Km '■= CIiZm-Mi- By Prohorov's theorem, each Mi is a relative compact subset of ty(Dp). Then we 
derive via Chebycheff's inequality for any Qn £ ^Pa^ 

N 

= e- 2M Eg„, N (exp(2(^l K « e ^})) 

< e~ 2Nl II + e 2i2 (e" 2 ' 2 (e ; - i)))^ (s ' w) < e "^. 

Therefore, 

oo oo 

pg jviiV (L i v ^m) < E £ M i) < E e " w ^ 2e_MAV2 ' ( 3 - 57 ) 

Z=M Z=M 

which implies the exponential tightness. 

□ 

4. Appendix 



4.1 A Lemma for pair measures 

In this appendix we provide a lemma for a unique characterisation of a pair measure with equal 
marginals via a vector in the coordinate space. Let E G N and v := E 2 — E in the following. The 
following lemma characterises pair probability measures Q £ • • • , E} 2 ) which have equal first 

and second marginal. Recall the definition for the first and second marginal, respectively = 
Ef=i Qti, k) and Q^(j) = £? =1 for any j G {1, . . . , E}. 

Lemma 4.1. Xei x = . . . ,x M ) G [0, l] 1 ' be a vector with the following properties 

E-2 

(1) x(kE+1) ^ 

fe=0 
(fe+l)E 

(2) < x (fcE+1) for k = 0, . . . , E - 2; 

(3) x (^+m < x (fc«+i) / or fc = o, . . . , E - 2; 

3=0 

E-2 E-2 



(4) 1 > (lx {kE+l) - x (fe+2+jB) ) . 



fc=0 3=0 

T/ien x defines uniquely a pair probability measure on {1, ... , E} 2 with equal first and second mar- 
ginal. 
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Proof. Let x G [0, l] u with the properties (l)-(4) in (4.1) be given. We define (E — l) 2 entries of a 
function Q: {1, . . . , E} 2 — > [0, 1] and E 1 entries of a function Q (1) : {1, . . . , E} — > [0, 1] in the following 
way: 

Q (1) (l):=i (1) ; Q(l,l) :=x< 2 >; ...; Q(l,£ l - 1) := x {E) 

Q^(2) :=x^; Q(2, 1) := s< B + 2 >; Q(2, £ - 1) := x< 2E > 

; (4.2) 



QW(E- 1) :=x« B - 1 > Q(£7 - 1, 1) := x« B -V ...; Q(£ - 1, E - 1) := x<"> 
respectively the £7 entries 

:=X (0-W) fo r j = 1,...,E-1 and Q (1) (£) := 1 - £ Q (1) (j). (4.3) 

j'=i 

The remaining (2.E — 1) entries of the function Q are functions of the given entries in (4.2) and (4.3) 
and are given as 

E-l E-l 

Q(j,E) := QW(j) - £ QU,k) and Q(E,j) := Q«(j) - £ Q(k,j) for j = 1, . . . , E - 1; (4.4) 
fc=i fc=i 

and 

Q(E,E) := 1 - £ fe). (4.5) 

0,fc)e{i,...,} 2 , 

The entries of the function Q defined (4.3) and (4.4) are elements in [0, 1] because of the properties 
(1), (2) and (3) of (4.1). Property (4) of (4.1) gives that Q(E, E) G [0, 1], because of 

1- £ QU,k) = l-{j22Q^(j)-J2Q(k,j)). 

(j,fc)e{i,...,} 2 , j=l j,k=l 

(x,y)jt(E,E) 

Thus Q(j, k) G [0, 1] for all j, k G {1, . . . ,E} and Q«(j) G [0, 1], j G {1, . . . , E} and £f=i QW(j) = 1. 
Clearly, 

E-l E-l E-l 

£ Q(j,k) = £ QU,k) + £ ( 2 Q W U) ~ £ {QU,k) + Q(k,j)))+Q(E,E) = 1. 
i,fce{i,...,£} fc,i=i j=i fc=i 

For any k = 1, . . . , E — lwe have from the second equality in (4.4) 

E-l 

QW(k) = Q(E, fc) + £ Qti, k) = Q w (k). 
3=1 

From (4.3), the first line of (4.4) and (4.5) we get 

E-l 

Q' 2 \E) = Q(E,E) + ^2Q(j,E) = 
3=1 

E-l E-l E-l E-l 

= !-(£ 2Q«(j) - £ J)) + £ (q (1) (j) - £ QV, k)) 

3=1 j,k=l 3=1 k=l 

= i-J2Q w (j) = Q w (E). 

3=1 
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Hence, altogether the defined function Q: {l,...,E} 2 — > [0,1] is a probability measure on 
{1, . . . , E} 2 with equal first and second marginal. □ 

4.2 Entropy estimation The following lemma gives an entropy estimation. As we are not aware of 
any reference for this we include the proof of this lemma. 

Lemma 4.2 (Entropy estimation). Let E be a finite set and P,Q € VP(E) two sub-probability 
measures such that 

^|P(x)-Q(x)|<a<l forae(0,±}. (4.6) 

x<=E 

Then the difference in the entropies is bounded as 

\H(P)-H(Q)\<-a\og^. (4.7) 

Proof. We write A(x) = \P(x) — Q(x)\,x € E. Since the function f(x) = —xlogx is concave and 
/(0) = f(l) = 0, we have for every < x < 1 — A, < A < ^ the estimation 

\f(x) - f(x + A)| < max{/(A), /(l - A)} = - A log A. 
Hence for < a < \ we get 

\H(P) - H(Q)\ <^\f(P(x)) - f(Q(x))\ 

x£E 

<alog \E\ — a log a, 

where we used for the last estimation the usual entropy estimate. □ 
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